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Introduction 


This handbook has been written to accompany the third-year module in 
Electromagnetic Engineering, given at the University of Cape Town. The 
main motivation for providing these notes in handbook form has been that 
students come to the course with a very wide variety of educational back¬ 
grounds; in particular, there are considerable differences in the amount of 
formal physics that they have covered. The notes for the course therefore 
aim to provide a description of electromagnetic field theory which can be 
used either as a first introduction or as a coherent revision programme for 
the subject. The emphasis is necessarily on the use of electromagnetic field 
theory by electrical engineers, and so the topics chosen for coverage have 
always been with that slant in mind. 

Throughout the notes, it is assumed that students have a working familiarity 
with the vector algebra that is studied as a core course in the second year 
at UCT. It is not assumed that students have taken second year physics, 
however. Only the first-year physics course is a prerequisite. 

The notes are presented in an open format, to encourage the writing of further 
comments in the margins. This is in the belief that it is better for students 
to be able to concentrate on the material during lectures, rather than feeling 
that they need to write copiously. The basic outlines of the lectures are laid 
down in the notes, and students are invited to read up in advance of each 
lecture, and to add to the handbook only where they feel it necessary. Much 
of this should preferably be done after the lecture, with reference to some of 
the many textbooks that are available on the subject. 

At the back of the handbook there are a few test and examination questions, 
which were given as part of the course in 1995 and 1996. They are included 
as a sample, indicating the standard required at the end of the course. 

I am very grateful to Richard Lord for helping to compile the text in RTj^X, 
and for producing the diagrams. The material for these notes is derived from 
several textbooks, but errors in interpretation are, of course, my own. I am 
always very grateful to be informed of any errors encountered in the use of 
the handbook so that they can be corrected in later editions. 
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To the Student 


This set of notes is designed to help you in your study of electromagnetic 
fields as part of the third year course at UCT. It is very important, however, 
that you appreciate the spirit in which the notes were produced: they are not 
intended to substitute for lectures, and they should stimulate you to look for 
more information in textbooks. Textbooks will give you much more informa¬ 
tion than is possible in a short handbook of this nature, and will enrich your 
understanding of the subject, particularly if you make a point of dipping into 
several works of references from time to time. 

It is also important for you to study this material on a regular and continuous 
basis, rather than to try to learn it in short bursts before tests! These notes 
are designed to help you to work steadily on improving your understanding, 
going forward at the rate of about one chapter per lecture. You may find 
that you need to brush up on your second-year maths. Put some time into 
this early in the course, if necessary—it will enhance your enjoyment of the 
later chapters! 

Field theory is a highly theoretical part of the background education of any 
electrical engineer. The fundamentals of field theory help us to understand 
many phenomena at a basic, physical level. However, there can be few prac¬ 
tical experiments in a course of this nature, and students often find it difficult 
to come to terms with the large amount of theory. Persevere, and the course 
will offer you explanations of many matters that are relevant to electrical 
engineers. Among these, we will tackle the following questions: 

• What exactly is meant by voltage and current? 

• What is the theoretical basis of capacitance and inductance, and how 
can they be calculated? 

• How do electric and magnetic fields store and release energy? 

• Do Maxwell's equations describe all that we know about electromag¬ 
netic phenomena? 
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• Why do electromagnetic waves travel through a vacuum at a velocity 
given by the well-known constant c? 

• Why are metals shiny, and why are radio waves reflected from layers in 
the ionosphere? 

As you can see from this short list, we will address some very fundamental 
issues in electrical engineering. I hope that you will enjoy using this handbook 
as part of your studies in this area, and that the skills which you gain in taking 
this course will be of value to you later in your career. 

G.M.T. 
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Electrostatics 







Chapter 1 
Electric Fields 


1.1 What is a field? 

• A wav to visualise a measured quantity mathematically, in a region 
with several dimensions. 

• e.g. Temperature in a room - a 3D scalar field; 
e.g. Pressure on the blackboard - a 2D scalar field; 
e.g. Velocity of a river - a 3D vector field. 

• Typically, some other physical quantity is the source of the field. Thus: 

- Sources of heat give rise to the temperature in the room; 

- Writing on the blackboard is a source of pressure; 

— Changes in the height of a river’s bed determine its velocity; 

— Mass is the source of gravitational fields in space. 

• The source of electric field is charge. 

1.2 The force between charges 


• Everyone has seen the force between charges in action: 

- A comb attracts small pieces of paper in dry conditions; 


9 


CHAPTER 1. ELECTRIC FIELDS 


10 


— A balloon sticking to a wall; 

— Xerography. 

• The fundamental experimental observation of electromagnetic theory: 

— Like charges repel and unlike charges attract; 

— The force between two charges in a vacuum was found by Coulomb 
to be 

f = (in magnitude) or f = K— |-r (in vector form) 

t 0 r z e Q r l 

(I- 

where 

* q\ ,q 2 are the magnitudes of the two charges in coulombs; 

* r is their distance apart in metres; 

* r is the unit vector from one charge to the other; 

* e Q (from experiment) is 8.854 x 10 -12 ~ ^ x 10 -9 F/m; 

* K is chosen to be -A to ensure that / is in newtons; 

* In the more general case where the charges are embedded 
in a medium, t is found from e 0 e r , where e r depends on the 
material separating the charges and is called the relative per¬ 
mittivity or dielectric constant of the material. 



1.3 Electric field 


• In the vicinity of a charge (or group of charges) we are often interested 
to know how other charged objects will behave, e.g. The behaviour of 
an electron in a TV tube with highly-charged X and Y plates. 


• To do this, we visualise the force per unit charge for each point in the 
region of interest: this is called the electric field. Thus: 


a 

Aq 


( 1 . 2 ) 


for a tiny test charge , A< 7 , too small to contribute significantly to the 
field in its own right. 
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• Hence, the electric field of a single point of charge, q, is given by 

NJc'' 


E = I< j$s. r = 


er 2 Aq 4r er 2 


\J $ v- 


_ | 


(1.3) 


Note: Since r is a unit radius vector from < 7 , the E-field points away 
from q in all directions, if q is a positive charge, or towards it if q is 
negative. 

• Recalling that charge, capacitance and voltage are related by Q = CV, 
we see from the previous two equations that the units of the electric 
field are either N/C or (C)(m/F)(l/m 2 ) =V/m. 


1.4 Superposition of fields 


• A principle of superposition holds for finding the electric field due to 
two or more charges - simply add the fields from each individual charge 
vectorially. 


• Thus, the field due to two point charges qi and q 2 is simply 

E = 


<72 . 

—J r i + - ~r 2 


47 rer( 1 ‘ \ieer\ 
where denotes vector addition. 


(1.4) 


• We often want to know the field due to whole lines or sheets or clouds 
of charge. To do this, we need to know the charge density , pi (in C/m) 
or ps (in C/m 2 ) or just p (in C/m 3 ) and the element of the group of 
charges, dl or dS or dV. We must then integrate the field due to each 
elemental charge, p\dl , psdS or pdV as appropriate. 
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e.g. The field of an infinite line charge 



Let the charge density all along the line be pi 


The elemental charge is therefore dQ = p\dz 


• For a test charge at the general point P, Coulomb’s Law gives: 


di dQAQ A _ dQ ^ _ dQ re r - ze z _ n M 
A Q 4ire c R 2 &Q 4ire 0 R 2 4ire 0 R 2 yV 2 + z 2 1 ' ’ 


• Integrating along the full length of the line of charge, all z-components 
cancel, so (using integration tables) the total E-field is calculated as: 





pirdz 

4ire 0 (r 2 + c 2 ) 3 / 2 r 


Pir 

4tT6 0 


z 

r 2 \fr 2 + ~ 2 


OO 


e r 


Pi . 

--e 

Z7 T€ a r 


— oc 


( 1 . 6 ) 
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1.5 Visualising an electric field 

• When you get answers like the above, it is important to interpret them. 
We can do this with a sketch of the E-field. 


e.g.l E-field of a single point : 


Airer 2 

Remember that the E-field is defined as the 
force per unit charge acting on a positive 
test charge. 




e.g.2 E-field of an infinite line charge : 

E =2^ 

Note that the field is not inverse-square. 
e.g.3 E-field of a pair or group of charges : 


Here you have to draw electric field lines 
that are everywhere tangential to the total 
electric field vector. 


* a 






\ r 


• Try experimenting with the iron-filings program to see how to arrive at 
the sketch of the field in example 3 above. Remember, however, that 
we are dealing with electrostatic fields, not magnetic fields, here - so 
the concept of iron filings is only an analogy. 


1.6 Fields of standard charge configurations 

Apart from the infinite line charge, there are other standard charge configu¬ 
rations whose fields you should be able to calculate and draw. You should do 
your own research (work in small groups) to make sure that you can clearly 
demonstrate the calculations leading to the results on the following page. 

• Can you see how 3 —> 2 as a —> oo? 
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Configuration 

E-field 

Sketch of the configuration 

1. A ring of charge of 


4 

* 

radius a, in the plane 

PldZ 

2e{a 2 + z 2 ) 3 / 2 * 2 



perpendicular to 



-^21 } 

the z-axis 

(on the z-axis) 



2 . An infinite plane 




of charge, filling 

Ps. 

-p 

x t 

S- 

2 e * 

X 1 

Z / 

the x-y plane 

■ y i 


3. A uniformly-charged 




circular disc of 

Psh ( , A A 

2 e \\/a 2 + h 2 hj 2 

4 

r*o i 

( r ,*, M 

radius a, in the plane 




perpendicular to 

(at (0,<j),h) on the z-axis) 



the z-axis 





^Ci . ,a 
Hi - 1 


















Chapter 2 

Electric Flux Density 

2.1 Electric Flux and Flux Density 

• So far we have dealt with E in V/m, which is strictly called the Electric 
Field Intensity. We have drawn field lines to illustrate electric fields. 
We need two more very important concepts for a much better way to 
handle fields. 

• Electric Flux ($) originates on positive charges and we visualise it as 
flowing along electric field lines towards negative charges (a sort of 
electric stream). 

Note: Flux is just another way of thinking about charge (1C of charge 
= 1C of flux or 'F = Q). If this seems strange, think of light, which we 
often conceptualise in the same way. 

• Electric flux is a scalar , but we are often interested in how many flux 
lines flow through a given area, and which way they are flowing. 

If near P in the small area dS , normal to a n , flux passes through dS 
parallel to a„, we then think of the electric flux density at P as being 

D = — an (in C/m 2 ) (2.1) 
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2.2 Gauss’s Law 


• Now consider an arbitrary charge Q. It may be a point charge or some 
linear, surface or volume distribution—it doesn’t matter. We surround 
it with some arbitrary surface, S, whose actual shape doesn’t matter 
either (we call this a gaussian surface). 

• At an arbitrary point P on 5, the flux density D through dS near P 
will almost certainly not be parallel to a n , the normal to dS. Instead 
a n and D will have some angle, 9, between them. 



( 2 . 2 ) 

= DtiSa n since D-a n = |D||an| COS 9 = |D1111 cos 9 (2.3) 
= D • dS (2.4) 

where dS is the vector surface element which has magnitude dS and is 
in the direction of a n . 

Integrating over the whole closed surface S now gives us 


V = fD-dS 


and, since = Q, we get 

^D-dS = Q 

or, for a charge distribution 

It>-dS = J pdV 


(2.5) 


( 2 . 6 ) 


(2.7) 


This remarkable result, known as Gauss's Law , says that the total flux 
out of a closed surface is equal to the charge enclosed by the surface. 
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2.3 Advantage of Gauss’s Law 

• Gauss’s Law enables us to calculate electric fields using the electric flux 
density , D, and not electric field intensity , E, directly. We look quickly 
at why this is an advantage. 

• Consider a charge Q at the centre of a sphere S (a gaussian surface of 
radius r). From Gauss’s Law, 



Q = 


fv.ds 

= D j dS 

= D(47rr 2 ) 


Q 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


4 xr 2 


where e r is the normal to dS on the 


sphere (and may also be written as the unit radial vector r). 

We also know that E = —^-e r from the previous lecture, so evidently 

4 Ter 2 

in this case, as in every other, 


D = eE 


( 2 . 11 ) 


where e is the permittivity of the medium. 

• In the vicinity of a charge, therefore, D is independent of the medium 
(or dielectric) but E is not. In problems involving multiple dielectrics 
it is a great advantage to calculate D from the given charge and then 
convert to E in the end dielectric. 

• Remember: Electric field , E, and Electric Flux Density, D, are both 
vector fields that have the same form and differ only by the constant 
factor e, determined by the medium. 
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2.4 Polarisation of the Medium 

• Before starting to use Gauss’s Law, we ought to understand why the 
electric field around a charge is affected by the dielectric medium, and 
how we quantify the effect. 

• Many media consist of atoms or molecules that are polar, they have 
distinct positive and negative regions (like the nucleus and electrons 
of an atom). A charge placed in such a medium tends to attract the 
regions of unlike charge, twisting or distorting the atoms or molecules 
of the medium. This is called polarisation. 

Note that polarisation tends to counteract the field due to the point 
charge. 

• We account for polarisation by writing D — e 0 E + P, where P measures 
the extra effect on D of this electric polarisation. Furthermore, in most 
materials the polarisation, P, is linearly related to the applied E-field. 
We choose to write this as: 

P = CoXeE (2.12) 

We then have the two equations: 

D = e 0 (l + X e)E (2.13) 

= eE (2.14) 


• Hence, since e = e r e 0 , we can write the relative permittivity of the 
medium as 


e r — 1 + Xe (2.15) 

The constant \t is called the electric susceptibility of the material of 
the medium. 
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e.g.l The field of an infinite line charge 


• Construct a special gaussian surface around the charge (i.e. a closed 
surface at which D from the charge is either normal and constant or 
tangential. Since the ^-component of D is zero in this case, the surface 
is a closed circular right cylinder: 

• By Gauss’s Law: 

Q= i D • dS = / D • dSi + / D • dS 2 + f D • dS 3 (2.16) 

JS J 1 J 2 J 3 

So 


piL = 0 + D{2nrL) + 0 (2.17) 


Therefore 


ID 


PiL 
2 tt rL 
Pj 

2irr 


(2.18) 

(2.19) 


Thus D = -^-e r and so E 
2 ?rr 



(as before but much more simply). 
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e.g.2 The field in a semiconductor depletion 

layer 


• An important and interesting practical example. Doped semiconductor 
in intimate contact with metal (or semiconductor of opposite doping 
as in the PN junction) gives rise to a depletion layer for 0 < x < d as 
shown. 


• In the x > d region, there is no net charge, and so Q — 0. Therefore 
E = 0 and hence D = 0 in this region, so all flux begins on the positive 
donor nuclei in the depletion layer and must end at the metal surface. 


• The special gaussian surface for this problem includes the plane x = x', 
where 0 < x' < d. The surface is closed, but only returns to x = x' 
at great distance from the depletion layer, so that effectively no flux 
passes through the curved return path of the surface. 

• By Gauss’s Law, Q = T> ■ dS so Npe(x — d) = D x (x). where Np 

represents the number of donor nuclei per unit volume in the depletion 
layer, and e is the charge on a single electron. Note that we have al¬ 
ready divided through by the cross-sectional area of the semiconductor 
device, and have now written the flux density as a function of x. 


Hence, the x-component of electric field is E x = 


Nne(x 

e 


d) 




CHAPTER 2. ELECTRIC FLUX DENSITY 


21 


e.g.3 The field between concentric spherical 

electrodes 

• Note that the dielectric medium between the spheres has permittivity 
ei for a < r < b, and permittivity for b < r < c. 

• The spherical symmetry implies that the electric field will be radial. 

• Let the charge on the inner sphere be Q and the charge on the outer 
sphere be — Q. 

• The special gaussian surface for the problem is a sphere, concentric 
with the other spheres, of radius a < r < b for medium 1 or b < r < c 
for medium 2. 

• By Gauss’s Law Q = <f D • dS = (DAiAirr 2 ). Hence D r = — - 

Js Airr 2 

anywhere between the two charged spheres. 

• The E-fields now follow easily: E r = —- for a < r < b and 

Airtir* 

Q 

E r = -- for b < r < c. So note that the D-field is continuous 

at r = 6, but that the E-field changes abruptly there. 





k 















Chapter 3 

Electrostatic Potential 

3.1 The conservative property 


• A charge, Q , placed in an electrostatic field, E, experiences a force 
F = QE. If the charge were to move within the field, it would therefore 
either gain energy or lose it, depending on whether it moved with or 
against the local value of the force. 

• In a field, F, whose strength varies in space, the differential work done 
on Q by the field as Q moves along vector dl is written 

dW = Fdl cosO = Fdl (3.1) 


Since total energy is conserved, this movement of Q also does work 
dW = —F • d\ on the system of charges that give rise to the field E. 


• Consider the work done in moving a test charge Aq from infinity to 
P, a point in the vicinity of a cluster of charges <?i, < 72,93 • • • at points 
Qi, QQ'6 ■ ■ •• Arising from the movement of A q, the total work done 
on the system (i.e. on q\,q 2 ,(} 3 - ■ •) is given on the next page, in terms 
of the diagram below. 
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W 


= z[-r 




/•pgi /-PQ2 y-pgs 

= - Fr-dl - F 2 d\ - F 3 -d\~... 

J oo «/oo . r «/oo 

y ' 


rC W) 


!o3 t>t f 6 ' 

,, <J ^ 


“ ~L 


PQ\ Aqqt cos d}dl 
47rerj 


-L 


A* v. 

P^2 cos d 2 dl 


47rero 




A/ 

PQs cos 6 3 dl 
Airerl 


rA^i_l PQl + rAgg 2 1 PQa + rA<?g 3 1 p ° 3 
L47rerJ 00 L47rer 2 J 00 L47rer 3 J 0O 

Ag?i Aqq 2 Aqqz 


AirePQi 4nePQ 2 4nePQ 3 


+ ... (since dr\ = d/cos etc.) 

(3.2) 


• So the total work done on the system depends only on the final position. 
P and not on the path taken by the charge, A q. Hence, if the charge 
moves around any closed path, no net work is done, so we can write 


F • d\ = 0 and so finally / e • <n = o 


( 3 . 3 ) 


This is known as the conservative property of an electrostatic field. 


3.2 The electrostatic potential 


• When a charge Q moves from point B to point A inside an electrostatic 
field E, the work done per unit charge is called the potential of point 
A with respect to point B. This is variously written 

$A ~ $b = V AB = W/Q = - \ A E • d\ (3.4) 

JB 


and the unit of the potential is thus the joule/coulomb or volt. 
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• Note: 

1 . Strictly, we only define potential difference. The “potential at B'\ 

= Vb, is actually arbitrary, and a point of zero potential may 
be placed anywhere for convenience. 

2 . For three points, A, B and C within E, since E is conservative we 
have V A b = Vac ~ Vbc- 

3. By our defining equation, Vab = — Vba- 

4. Depending upon the co-ordinate system in use in a particular 
problem, the expressions for d\ are: 

(a) d\ = dxe x + dye y + dze z in cartesian co-ordinates; 

(b) d\ = dre r + rdcfe^ + dze z in cylindrical co-ordinates; 

(c) d\ = dre r + rd9eg + r sin Odcfe# in spherical co-ordinates. 

• If the potential at infinity is taken as zero (often but not always con¬ 
venient) then the potential of the test charge A <7 at the general point 
P is simply 

W __ , Si | Si q A 

A q AirePQi A’KtPQ 2 ^-KtPQ 2 ~ AireRi 

where PQ. t = Ri = |r - rj| = y/(x- x ') 2 + (y - y ') 2 + (z - z ') 2 
Note how the prime superscript denotes a source of electrostatic field 
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• If the source of the field is a continuous distribution of charge, of den¬ 
sity p{r ) within a volume V, the potential function of the distribution 
is 


*(r)= / 

h 


$(r) = f 


p(r')dV 
v AiteR 

p(r')dV‘ 


lv AireR 
(C is just a constant) 


if the potential at infinity is zero, or 
A C if zero potential is defined as being elsewhere. 


e.g.l The potential of a point charge, Q 


Since the field is known to be radially symmetric (E = ;jd2_r), we work in 
spherical co-ordinates with d\ replaced by using 


d\ = dr e r + r dO eg + r sin 6 d<f> 


The potential is then 


V. 


AB 


f A rr A rr A 

— I E • d\ = — / E r ■ e r dr = — 

J B J r q J rg 


E r dr 


rr/ 
Jr B 


r A Q 


Alter 2 


dr — 


[11 

rA 

Q 

[1 

1 • 

Lr. 

rB 

4ire 

ir A 

r B ■ 


(3.6) 


(3.7) 

(3.8) 


Then, if the reference point B is at infinity, we have the potential at any 
point A with radius vector r: 


*(r) 


Q 

Liter 


(3.9) 
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e.g.2 Potential near a line charge of density pi 


We know the field to be E = cylindrical co-ordinates, so dl is 

replaced by using 


dl = dr e r + r d<t> e $ + dz e z 

The potential is then 


(3,10) 


Vab 


[A rr ,i rr A 

/ E • dl = — / E T ■ e r dr = — E r dr 
J B Jtq J rg 


(3.11) 


IB 

£ 2 £r dr = ik |ril " = kk lln |rsl - ln (3 - 12) 


With the reference of zero potential at r B = oo we would get a nonsensical 
result, so we arbitrarily define zero potential as being at the fixed point r Q . 
Hence, $(r 0 ) = 0 which gives us: 



V° B = ^ [In |r B | - In |r 0 |] = 0 

(3.13) 

so 



and so, in general: 

ln|r B | = In |r 0 | 

(3.14) 

*(r) = 

27re t ln l r l ln l r °l] = 27re ^ ^ ° 1 

(3.15) 


where, again, C is just a constant. 
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e.g.3 Potential of an infinite sheet of density p$ 


We first need an expression for the electric field near the sheet. One way 
to get this is to construct the cylindrical special gaussian surface shown in 
the diagram. Then, by Gauss’s Law, the small enclosed charge q is given by 
/ 5 D ■ dS — q. We can imagine the length of the cylinder to be arbitrarily 
short, so that flux flows through both flat ends of the gaussian cylinder. 
Thus: 

q = Dj • dSi + J D 2 • dS 2 = Da + Da (3.16) 

Hence, D = q/2a. 


Since the position chosen for the gaussian cylinder was in no way special, 
and since q/a — ps, we can write for the whole sheet that D = ps /2 and so 


E = ±ge, (3.17) 

with the E-field acting only in the indirections. Then, setting the reference 
potential at the sheet (at z — 0) as $(zb) = K, we can find the potential: 


V AB = - [ A E • dl = - \ ZA E z • e z dz = 

J B J zb 

at a perpendicular distance za from the sheet of charge in the positive z- 
direction. 


Z B 


T T PS 
= - 2^ 


(3.18) 














CHAPTER 3. ELECTROSTATIC POTENTIAL 


28 


3.3 The gradient of the potential 

• In the vicinity of charges, we now know what is meant by the electric 
field and by the potential function. If you join points of equal potential 
in an electrostatic field, you get what are known as equipotential lines. 
We now demonstrate that these lines of equipotential are always at 
right angles to the field lines. 

• Since $ = — /E • dl. d$ = —E • dl. In cartesian co-ordinates, we can 
replace dl with xdx + y dy + zdz and E with E x Se + E y y + E z z (note 
the slight change in notation), so 


d$ = —E x dx — Eydy — E z dz (3.19) 

Now too, is just a function of x, y and z, so we may write (see 
MAM280W): 


,, 5$ . d$ , <9$ , 

d$ = —da; + -z-dy + — dz 
ox dy dz 


(3.20) 


from which, clearly 


E x = 


dx 


Ey 


_ _3<f> 

d y tz ~ dz 


(3.21) 


Hence, we can write the field E in terms of the potential only by stating 

-V$ (3.22) 

where the vector operator V should be very familiar to you! 


E = -hte +y w +z & 


• Now, since V$ is a vector showing the magnitude and direction of 
the greatest variation in potential seen from a point, the field lines lie 
in that direction too (E = — V$). Equipotentials are the lines of least 
variation in potential, along which d$ = V$-dl = 0, and this condition 
will only apply if V3> L dl. Hence, the equipotentials are everywhere 
perpendicular to the field lines, which was what we set out to prove. 
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3.4 Standard charge configuration potentials 

The table below gives some standard charge configurations of high symmetry 
for which you might like to calculate the potential functions. You will prob¬ 
ably need to research the methods for doing these calculations, and should 
work in groups and look for texts that give guidance. Make sure that you 
are each able to perform these derivations independently, having done the 
research together. 


Find the potential of: 

Sketch of the configuration 

1. A sphere of charge Q 

in a vacuum, at a distance 
r from the centre of the sphere 
where r is greater than 
the sphere’s radius 

0 

2. A circular disc of charge 
of radius R and charge density 
ps , at a distance z along 
the z-axis which is the 
perpendicular axis of the disc 

J 


/ 



3. A circular ring of charge 
of radius R and charge density 
pu at a distance z along 
the z-axis which is the 
perpendicular axis of the ring 

t 

> 

c* 11 


y 


















Chapter 4 

The Field Between Electrodes 


4.1 Graphical field mapping 

• We now develop a technique for drawing the field and equipotential lines 
in the vicinity of charged electrodes of arbitrary shape. Our drawings 
will, of course, be 2-dimensional, but they are a good approximation 
for more realistic cases. 

• We start with a known potential difference between two electrodes 
drawn in a diagram. Divide this difference into a convenient num¬ 
ber of potential levels, and add to the diagram (using commonsense 
and some guesswork) an equipotential line for each level. 

• It is best to start in regions where the field is straightforward and 
uniform. Thereafter, note that the equipotentials will tend to hug 
sharp, acute angles in the electrode boundaries and to avoid regions of 
reflex boundary curvature. 

• Then draw in the field lines, perhaps using a different colour. Remem¬ 
ber that the field lines are everywhere orthogonal to the equipotentials. 
Ensure that the resulting “curvilinear rectangles” are all similar in their 
side-size ratio. This will mean that you are drawing regularly-spaced 
lines of potential and field in your plot. 
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The proof of the above statement is as follows. 

In the diagram of a single curvilinear rectangle in a 
plotted field, An is the distance between two 
adjacent equipotentials, and As separates two 
adjacent field lines. The field strength in this small 
region is therefore AO/An V/m. Now the flux 
flowing between the field lines can be written as 

eAOAs 

AO = DAs = eEAs = —- 

An 



(note that we are so far ignoring the depth of the field into the paper). 


From this we have 


As 

An 


AO 

eAO 


(4.2) 


so if the flux per tube, AO, the potential difference per division, AO, 
and the permittivity, e, are all constant, then the side ratio, As/An 
will be constant too. 


• If your rectangles are not all similar, then your first guesses about the 
positions of the equipotentials were at fault. Do not be afraid to re¬ 
draw the diagram several times, until you are satisfied with the result. 


4.2 Field strength 

• You can read the field strength straight from a field map. For example, 
if V 0 in the drawing on the next page is 100 kV and d = 100 m, then 
E — 10 5 /10 2 = 1000 V/m at each edge of the drawing. Above the 
spike, however, this value is much greater because of the tendency of 
the equipotentials to hug the sharp angle. 
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e.g. Draw equipotentials and electric field lines 
between these objects 
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4.3 Capacitance 


• The capacitance between two electrodes of excess charge ±Q with a 
potential difference between them of <&a — is defined as 

Q 


c = 




(4.3) 


It is easy to sketch the equipotentials and 
field lines for a parallel plate capacitor. 
Ignoring for a moment the fringing flux that 
curves around the edge of the plates, we can 


1 




fm 


♦« 



*. !!!!!! 


ft A 

obtain, from Gauss’s Law, D = ps/2 for each plate. So with both plates 
considered, D = ps, the surface charge density. Then by definition of 
potential, and with the plates perpendicular to a cartesian z-axis: 




( A f A 

- = - / E • d\ = — / 

Jb Jb 


A D, 


zdz = 


Psz 


t 




Z B 


Psd 


(4.4) 


• If each plate is of area A then evidently Q = psA, so we get the well- 
known equation for parallel-plate capacitance: 

Q _ psA _ eA 

$a ~ $b Psd/e d 


C = 


(4.5) 


• Capacitance per unit length can be calculated from graphical field maps 
with two electrodes. The charge between the electrodes is the same as 
the flux ('L = C), and both are equal to the number of “flux tubes”, Nj, 
multiplied by the flux per tube, A'L. Similarly, the potential difference 
between the electrodes is just the number of potential divisions, N p , 
times the potential difference per division, A$. Hence 


C = 


Q 


$A - $ 


B 


Nf&<S> N f ft As 


7V P A$ 


N v \ An 


< 4 - 6 » 


assuming that the graphical field map consists of small squares so that 
As « An. Thus in the parallel-plate capacitor shown above, there are 
4 potential divisions and 26 flux tubes, so the capacitance is 


C ~ 10- 9 /36tt x 26/4 = 57.5pF/m 
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4.4 Conductance 

• Up to now, we have assumed that the electrodes in our plots were 
surrounded by a vacuum (e = t 0 ) or by an insulating dielectric (e = 
e 0 e T ). If, however, the surrounding medium is conductive , then the 
electric field E due to the charge on the electrodes will apply a force 
on the free electrons in the conducting medium. This will cause the 
electron's to move, and hence to establish a current. The current density 
(J in Amps/m 2 ) that results is related to the electric field by 


J = ctE 

where a is the conductivity of the medium. 
• It should be clear, considering a 
small region of the medium, that 


G> 





Af» # 


(4.7) 


-1 


so that 


A J — J As = crEAs = 
AJ 


crA3>As 

An 


ir&i 


As 

An 


<tA$ 


(4.8) 

(4.9) 


• The conductance per unit length between the electrodes is defined as 

G = 7 „ (in 0/m) (4.10) 

where the total current , /, is simply NfAJ (because the current flows 
along the flux lines). Hence 

r _ j _N f AJ N f Nj 

<$> A - N P A$ N p V AnA$ J ~ N p 1 ‘ ’ 

so long as the rectangles on the graphical field map are small squares 
for which As/An rs 1 


• We arrive finally at the useful relation 

G=-C (4.12) 

which you will meet in your work on transmission lines. Note that the 
rtsistanct of a conductive medium is just R = l/G (in 0). This gives 
us a way to find the resistance in a field map if a is known. 




















Chapter 5 


Divergence and Laplace’s 
Equation 

5.1 Integral and differential forms 


• We have seen electric field intensity and electric potential related by 
equations in both integral and differential form. Thus 

$ = - J E ■ dl and E = —V$ (5.1) 

• Gauss’s Law in integral form is f s D • dS = Q , where we recall that the 
charge, Q , may be written as Q = f v pdV. Differentiating both sides 
amounts to dividing through by AV and taking the limit as AV —* 0: 


Urn . lim kg. 

AV AV 


(5,2) 


We recognise the RHS of the above as simply p, and we propose to 
write the LHS as V • D. pronounced div D. Hence: 

VD = p (5.3) 

We call this the differential form of Gauss’s Law. 
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5.2 The meaning of divergence 

• Divergence refers to the net amount of a vector quantity per unit volume 
that emerges from an infinitesimal volume surrounding some point. It 
is therefore an expression of whether the quantity has a source at the 
point of interest. We can see what this means by considering a small 
cubic volume, AV , that lies within a field as shown. 




If the small cube lies inside an electrostatic field, 
it should be clear that the flux will be the same 
through any two of its opposite walls unless the 
flux density varies between two other walls, perp¬ 
endicular to the first two walls. In fact, for a 
small volume, the difference in any vector between 
two opposite faces is approximately the gradient of 
the function times the distance between faces. This 
approximation becomes exact in the limit as the 
sides of our small cube tend to zero length. 




• Let the small cube be centred at ( x,y,z ) and have dimensions Ax x 
Ay x Az. Considering the x-component of D, written D x , the above 
theorem gives us expressions for the flux density in the x-direction at 
the sides of the cube: 


D x 

D x 



D -(*) + f ^ 

Ax dD x (x) 
x[ ] 2 dx 


(5.4) 


• The flux out of the right-hand face, from the first flux density, is: 

Ax dD x (x)\ 


AyAzD x (x + ^) = AyAz (d x {x) + 


2 dx 


(5.5) 


and the flux into the left hand face is: 

AyAzD x (x - = AyAz ^D x (x) 


Ax dD x (x)\ 
~2 dx J 


5 (r) 


(5.6) 
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Hence the net flux out of the cubic volume in the x-direction is the 
difference of these two: 


Net flux out (x direction) = AyAzAx 


dDx 

dx 


(5.7) 


You can then argue similarly for the y- and z- components of D to see 
that the total flux out of the small volume is: 

d D d /~) c) J~) 

Total flux out = AxAuAz - - + AxAyAz - - + AxAyAz — - (5.8) 

ox uy oz 

Now, by Gauss's Law, the total flux out of our small volume must equal 
the charge that it encloses, which is just pAxAyAz, so in the limit we 
can now write: 

^ + = „ (5.9) 

ox oy oz 


But we can define div D as the vector operation: 

div D = V • D = 


- 


v & / D - 

and so we again conclude that V • D = p 


Dx | _ dD x dD y 

Uy 1 “ dx 


+ 


dy 


+ 


dD z 

dz 


(5.10) 


• Note the compactness of this elegant expression: it says that “the 
source of flux density is charge”, something we began with back in 
lecture 1, but which we now see to be another aspect of Gauss’s Law. 
Note also that you can get back to the integral form of Gauss’s Law by 
using the divergence theorem , a mathematical identity that holds for 
any vector. F: 

-FdV = dS (5.11) 


Integrating p = V • D and using this theorem, we get 

J pdV = fv T)dV = j D dS = Q 


(5.12) 
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• YD can also be defined in cylindrical and in spherical co-ordinates, 
using arguments much like the one presented above for the cartesian 
case. You should make sure that you have researched these arguments, 
and are capable of showing the following from first principles: 


V D 

YD 


hl-y D ’ )+ 


i d 

r sin 9 89 


(Do sin 9) + 


1 8D 




r sin 9 8<t> 


Id 1 dDb 8D Z 

~a~( rDr + + ~!T~ 

r dr r d(j) dz 


(sph.) 

(5.13) 


e.g.l Find V • D for a solid sphere of charge 
density p and radius a, embedded in a vacuum 


Inside the sphere, using a gaussian sphere of radius r < a, Gauss’s Law gives 

<f D • dS = Q so DAnr 2 = -7r r 3 p so D = — so D = — e r (5.14) 
Js 3 3 3 

Then, using spherical co-ordinates, 

V.D-14-fr*^-!^-, (5.15) 


" 2 dr 

exactly as expected in a region of charge density p. 

Outside the sphere we work with a gaussian sphere of radius r > a to get 


(f D • dS = Q so DAnr 2 = -7r a 3 p so D — so D = —ye r 

Js 3 3 r l or 1 


Then, using spherical co-ordinates. 


V D = 4— (r 2 —) = 0 
r 2 dr l 3r‘ / 


(5.16) 


(5.17) 


again as expected in a region of zero charge density. 
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5.3 The equations of Poisson and Laplace 


• Using the three relations V • D = p, E = — V5> and D = eE, we can 
now T write 

V • E = - = —V • V$ (5.18) 

t 

and this is generally seen as a differential equation called Poisson’s 
equation , which is given below in several forms: 


V • = V 2 $ 


d 2 $ d 2 <$ d 2 $ 

'd? + ~df + ~d? 


(5.19) 


where V 2 $ = div(grad $) is called the Laplacian of $. 

• In regions of no charge, another important equation, called Laplace’s 
equation , results: 

V 2 $ = 0 (5.20) 

Solving Laplace’s equation to find the potential function $ amounts to 
solving a second order, partial differential equation in which $ may vary 
in 1, 2 or 3 dimensions, and may be expressed in cartesian, cylindrical or 
spherical co-ordinates. As you can imagine, the mathematics can vary 
from very easy through to rather complicated! We restrict ourselves to 
some very simple examples below. 


• In order that you can be comfortable with Laplace’s equation in cylin¬ 
drical and spherical co-ordinates, you should, in your own time, take 
the expressions that you derived for V • D in those co-ordinate systems 
and show that Laplace’s equation becomes 


V 2 $ = 


V 2 $ = 



2 d$\ 1 d ( . -d*\ 1 d 2 $ 

r ~ ' + ° • ede[ sm0 del + r 2 sm 2 ed6 2 


dr 


r 2 sin 


i- ( r 2*) + 

r dr \ dr ) r 2 dcf) 2 


1 d 2 $ d 2 $ 

+ dz 2 


MO 


(sph.) 

(5.21) 
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e.g.2 Potential between parallel conducting 

plates 


If the region between the parallel plates 
is charge-free, then we can write: 


d 2 $ { d 2 <S> ( d 2 $ 
~d^ + ~d^ + ~d^ 



Neglecting fringing, the potential is a function of z only , so = 0, from 
which we can integrate once to get ^ = A and integrate again to obtain 


$ = Az + B 


(5.23) 


We can evaluate the constants A and B by using the known values of the 
potential given to us as the voltages of the plates. Using this information is 
known as applying the boundary conditions. 


So, since at z = 0 we know $ = 0, we can see that B = 0 

and since at 2 = d we know $ = 100, we also obtain A = 100/d. Thus 

$ = 100 (volts) (5.24) 


This gives us an expression for the potential at any point (x,y,z) between 
the plates, and it suggests that equipotential levels would be equally-spaced, 
which we know to be the case. We can now find the expression for the electric 
field intensity, like this: 


E = -V$ = - 



_ _d_ 

dy^ dz 7, dz 




(5.25) 


and you should recognise this as having exactly the magnitude and direction 
that we would expect for an ideal, parallel-plate capacitor 
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e.g.3 Potential between coaxial cylinders 


If the region between the coaxial 
cylinders is charge-free, then we 
can write Laplace’s equation in 
cylindrical co-ordinates: 



\_d_ f 1 d 2 <L d 2 <S> 

r dr y dr ) + r 2 d<t> 2 + dz 2 


(5.26) 


Assuming no variation in potential with changing angle <£, and also assuming 
that the cylinder is long enough that there is no variation of potential with 
2 , Laplace’s equation reduces to 


i d ( d9\ n 
r dr \ dr J ~ 


so 



d$ 

so r—— = A 
dr 


so —— — Air so $ = Ain |r| + i? (5.27) 

dr 

Next we apply the boundary conditions: 

Since at r = a we know $ = 0, we can see that A In |a| + B = 0 and 
since at r = b we know $ = V 0 , we also get A In |6| + B — V 0 , 
and we can then easily obtain 


A = 


Vo 

In |6| — In |a 


and 


— Vo In jaj 

(In \b\ — In |a|) 


(5.28) 


The expression $ = A In |r| + B can now be used to find the potential at any 
radius r. knowing only the dimensions a and b of the coaxial cylinders, and 
the voltage on the outer cylinder with respect to the inner. 




Chapter 6 

The Uniqueness Theorem and 
Energy 

6.1 The uniqueness theorem 


• There is an important theorem which states that if we are given a 
region where Laplace’s equation holds, and if we know the boundary 
conditions, then the solution to the equation is the unique expression 
of the potential of the field in that region. 

• We can prove this uniqueness theorem by contradiction. Imagine that 
we have a charge-free region, with potentials specified on its boundary. 
Let us also imagine that there are two solutions for this potential, 
written $1 and $ 2 - At the boundary, both solutions must give the 
same potential (the potential specified), and so we can write for the 
boundary: 


= $ 2 SO - $ 2 = 0 


( 6 . 1 ) 


• Both potential functions would have to be solutions of Laplace’s equa¬ 
tion, and so we can write for the whole region: 


= 0 and V 2 $ 2 = 0 and hence V 2 ($j - $ 2 ) = 0 (6.2) 
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• Now recall the divergence theorem : for any vector F 

J v V • FdV = j> F • dS (6.3) 

We choose to use the divergence theorem with F = ($i — $ 2 )V($i — $ 2 ). 
This gives 

J V • [($, - $ 2 )V($j - <S> 2 )]dV = ^[($1 - $ 2 )V($, - $ 2 )] • dS (6.4) 


• Another mathematical identity (which you may find it very instructive 
to prove to yourself) states that for a vector A and a scalar if, 




V • OJA) = »V • A + A • V» 


(6.5) 


^ f 1 elf X 


Noting that — <I> 2 is a scalar and V($i — $ 2 ) is a vector, the previous 
equation may now be written 


h \J . A ' ' 

dt j ( pfi ■ ' — / v (« >1 -^2)V 2 (§ 1 -$ 2 )dV+^[V($ 1 -$ 2 )] 2 dy = jf($ 1 -$ 2 )V($ 1 -$ 2 )-dS 

dz V V 5 ( 6 . 6 ) 




7>* 


\. V4 r 


• Now the first term of equation 6.6 refers to the whole region and equals 
p N4 zero (by equation 6.2), and the term on the RHS refers to the boundary 
J X surface, and equals zero (by equation 6.1). We are left with 


J [V($, - $ 2 )} 2 dV = 0 (6.7) 

The gradient of a real scalar being real, and therefore having a square 
which is positive or zero, the gradient V($i — $ 2 ) can only be zero. In¬ 
tegrating, we obtain a relationship between the two potential functions 
which must apply everywhere: 


$1 — $2 = constant 


( 6 . 8 ) 
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• But if this applies everywhere, then it applies at the boundary, where 
we know the two potentials are equal. We are obliged to conclude that 
the constant in the equation above is zero, and that the two potential 
functions are equal. Our initial assumption of two different potential 
solutions has led to a contradiction, and we have proved that solutions 
to Laplace's equation with boundary conditions given are unique. 

6.2 Energy in fields 

• Evidently, electric fields store energy (they are capable of accelerating 
charge q through a potential drop v, and the energy required to do this 
is w = qv). We now aim to find expressions for the electrostatic energy 
stored in a field, in terms of field quantities such as E, $ and D. 

• The energy required to assemble a group of n charges, bringing each 
charge from zero potential at infinite distance to a particular potential, 
<fq, in an arrangement or system of charges is 

C (6-9) 

where the factor of 1/2 compensates for a formula that would otherwise 
count the contribution of each charge twice : once as it is brought up to 
the other charges, and once as each of them is brought towards it. 

• If the charge distribution is continuous, with charge density p, we write 

U = \tp$dV (6.10) 

Z J V 

where the potential, $, is a function of the position within the volume 
over which the integral is performed. Recalling that p = V • D, this 
becomes 


u = 


( 6 . 11 ) 
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and then, since for any scalar and vector A we have 


V • (u’A) = ?/>V • A -f- A • so ip {V • A) = V • (i/>A) — A • Y/ip (6.12) 
(you have recently checked the proof of this for yourself), we can write 

f/ = I^V.($D)dF-i^D-(V$)dF (6.13) 


• The first term of this, by the divergence theorem, now reduces to 


!>• i ' , f 

1 / V- ($D)dV = l l $D-dS =0 (6.14) 

2 Jv 2 Js 


since we are interested in the total energy due to every part of all of 
the fields of the charges present, and this involves integrating over a 
surface at infinity, at which distance the product of D and $ will have 
shrunk to nothing compared with the area of the surface of integration. 


' 2 ^ 1 ®' 


• We are therefore left with 

U = -]- f D • (V$)dV = \ j D • EdV (6.15) 

Z i/V £ JV 

This equation confirms that total energy resides in the field represented 
by D or E, and we can use it to express the energy present in each small 
volume element of the region by differentiating to get 

dU = -D-E dV 
2 


(6.16) 
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• Considering now a parallel-plate capacitor, in which 

— the plates have area A and separation d ; 

— the voltage drop between the plates is V c ; 

— the electric field has constant value E = VJd\ 

— flux density D = tE = eV c /d is constant, as e is constant: 

the energy stored in the capacitor’s field (neglecting fringing) is 


U r . = 


-DE [ dV 

2 Jv 

5(f) (!)<-> 




(6.17) 

(6.18) 
(6.19) 


But we recognise from equation 4.5 that cA/d = C, the capacitance 
of the parallel-plate system; and so we have, just from field theory 
considerations, the well-known result 


Uc = ^CV c 2 


( 6 . 20 ) 






Chapter 7 

Electrostatic Boundary 
Conditions 


7.1 Electric field due to dc current 


• When a dc potential is applied to a conductor, a steady dc current will 
flow. The electrons that participate in the current give rise to an electric 
field which, although not strictly static, is of interest to us now because 
it does not vary in time. Such fields are called quasistatic fields, and 
analysis of them proceeds similarly to our study of electrostatic fields. 

• In an ohmic conductor (a conductor obeying Ohm’s Law), the current 
density , J, is related to the applied field E via the conductivity , a. Fur¬ 
thermore, the relationship between electric field and potential allows 
us to write 

J = cxE = (7.1) 

• A stationary (i.e. non-time-varying) current will remove exactly as 
much charge from any closed region as it adds to it, per unit of time. 
There can be no build-up of charge either spontaneously, or due to a 
stationary current. This principle of continuity can be written 

jfj-dS = 0 (7.2) 

or. dividing through by AV and taking the limit as AV —> 0, 

VJ = 0 (7.3) 
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• We therefore have V • (—<rV3>) = 0, so with a a constant we can write 

V • V4> = V 2 $ = 0 (7.4) 

In other words, the potential that is due to a stationary current satisfies 
Laplace’s equation, just like the potential due to static charge. 

• We have already seen that solving Laplace’s equation amounts to solv¬ 
ing a differential equation in several variables. The rest of this chapter 
investigates the boundary conditions that we use to evaluate the arbi¬ 
trary constants in such calculations. 


7.2 The boundary between dielectrics 

• Since the field that we wish to study is often confined to a region 
occupied by one material and surrounded by other, different materials, 
we often gain useful information from knowing how E, D, $ and J will 
behave at the boundary. 


• Consider an arbitrary boundary between 
dielectrics 1 and 2 of permittivity ej and 
€■ 2 , with a surface charge ps on the boundary. 
We already know that the normal flux density, 
D n , is the same in both materials if there 
is no surface charge. Thus, if ps = 0, the 
component of D normal to the boundary is 
equal on both sides of it, and we can write 

D n \ = E) n 2 



This is not the case, however, if a surface charge is present. 

• Imagine a cylindrical gaussian “pillbox” embedded in the two materials 
at the boundary as shown above. Let the area of each flat circular 
face of the cylinder be AS, and let the cylinder be so short that the 
tangential component of flux density passing through it is negligible. 
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The charge enclosed by the cylinder is psAS, so clearly Gauss’s Law 
gives 

D n iAS — D n2 AS = psAS (7.6) 

from which 

D n i — D n 2 = ps (7-7) 

• If ps = 0, then we can use D = eE to see that 

t\E n \ = t2E n 2 (7.8) 

and that therefore the normal component of electric field intensity 
changes at the boundary because of the dielectric discontinuity. Now 
imagine traversing the closed path that crosses the dielectric boundary 
in the diagram. Let the long sides of the rectangular path have length 
A/, and let the short sides be so short that their contributions to / E-dl 
as we traverse the path are negligible. This means in effect that the 
potentials on either side of the boundary are equal (i.e. = <f> 2 ). We 

know from the conservative property of the field that 

£ E • d\ — Et\Al — E n Al — 0 (7.9) 

from which 

E n = E a so — = — (7.10) 

Ci e 2 

• This vector diagram summarises 
the effect on E of a change 
in dielectric at the boundary. 


It should be clear that 

9i = tan -1 (-zp-) an d #2 = tan -1 (^~) (7-11) 

V E n i / \ E n 2 / 

so, by equations 7.8 and 7.10, 

02 — tan -1 ( 77 ——^ = tan -1 (— tan^i^ (7.12) 

\E n i Ci / Vej ) 

and so if e 2 > ej, then the total E-vector is deflected away from the 
normal at the boundary. 
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7.3 The boundary between conductors 

• At the boundary between conductors, the normal component of cur¬ 
rent density, J n , is unchanged across the boundary: this is a necessary 
continuity, because there would otherwise be a build-up of charge at 
the boundary. Hence 

(J 2 

Jn 1 = Jn 2 SO CT\E n ] = <X 2 .E n2 &nd En\ — E n 2 (7-13) 

<?1 




The tangential electric field is unchanged 
across the boundary for exactly the reason 
that we saw when considering two dielectrics. 
Hence E n = E t 2 - Again, the angles and 62 
can be determined from the conductivities of 
the materials at the boundary, and equations 
very similar to equations 7.11 and 7.12 will 
be obtained. In particular, 

, / (J2 ^ 

O 2 = tan 1 ( — tan 

\<7l J 



Note that if material 1 were a semiconductor and material 2 a good 
conductor such as a metal, then C\ « cr 2 , so ^ —> 00 , so 0 2 —> 90°. 
We generally assume that metals are perfect conductors, in which both 
components of the electric field are zero. 

• At the boundary between a conductor and an insulator (a = 0), the 
normal component of current density, J n , must be zero on both sides 
of the boundary. Between a perfect conductor and an insulator, the 
tangential component of electric field, E t , must be zero on both sides 
of the boundary. For this reason, field lines always leave or meet a 
perfect conductor at right angles. Within a perfect conductor there can 
be no gradient of potential and therefore no electric field. 
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e.g. Example to illustrate boundary condition 

rules 


List the boundary conditions and draw the field lines for the structure below. 
You can assume that cr 2 >> G\ and ti > e Q , with <r 0 = 0. 
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Chapter 8 

Magnetostatic Fields 

8.1 Introduction 


• Magnets (lodestone) have been known for at least 800 years, and their 
use has had a profound effect on human history. 

• Just as electric charges are either “positive” or “negative”, magnetic 
“charges” or poles are of two types, called north and south. Note that 
magnetic poles are always only found in north-south pairs, never singly. 

• Several basic observations about the behaviour of magnets have been 
made, and must be included in a full theory of how magnets work: 

- Magnetic material attracts iron (known from antiquity); 

- Magnets have poles of two types: like poles repel each other and 
unlike poles attract each other (from about 1600); 

- A magnetised compass needle is deflected by a nearby electric 
current (A heuristic discovery by Oersted in 1820); 

- A constant current can give rise to an electromagnet with very 
similar properties to a permanent magnet (also 1820); 

- Magnets will deflect an electric current flowing in a vacuum or 
semiconductor by the vacuum tube or Hall effect (20th century). 

• We can explain these phenomena with the concept of magnetic fields. 
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8.2 Static magnetic field 

• A permanent magnet or a steady current is associated with a magnetic 
field. This is justified because the magnet or electromagnet is seen to 
exert a force on a nearby current-carrying wire, for instance. We know 
from electrostatics that forces in several dimensions are well-handled 
by fields. 

• Whatever the source of the field (permanent magnet or electromagnet), 
we define a vector B, called the magnetic flux density , so that the force 
exerted by the field on a small current element of length dl and carrying 
current I is of magnitude 


df = Jjdl||B| sin 9 

where 6 is the angle between d\ and B. 



( 8 . 1 ) 


• Now df is at right angles to both dl and B, and determined by the 
advance of a right-hand screw when dl is rotated onto B through the 
smaller angle between them. This is conveniently written as 


df = 1 dl x B (8.2) 

(Recall that if C = A x B, then C is a vector of magnitude |A||B| sin 9 , 
where 9 is the smaller angle between A and B. The direction of C is 
determined by the advance of a right-hand screw rotating A onto B 
through 9.) 

• The magnetic flux density, B, is measured in tesla (T) or webtrs per 
square metre (Wbm -2 ) or even Vsm -2 . We also have use for another 
field quantity called the magnetic field intensity, H, which is related to 
B via a property of the medium called the permeability, p. Thus 


B = //H (8.3) 

an d H is measured in Am 1 while p is in henries per metre, Hm h 
The permeability of free space is p a = 47r x 10 -7 Hm -1 . 
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8.3 


The Biot-Savart Law 



• Immediately following Oersted’s discovery of the magnetic effect of a 
steady current. Biot. Savart and Ampere deduced the relationship that 
expresses the size and direction of the magnetic field intensity. 

• The effect of a current / flowing through a small line element of length 
d\ at a point in space defined as r with respect to an arbitrary origin 
can be stated at any other point, P. If P is at the point defined by r 
from the origin, then the magnetic field intensity at P due to the flow 
of I in d\' is given in magnitude by: 


dH{ r) 


l\r')dX sin </> 
4ttR 2 


(8.4) 


where R = |r — r | and <j) is the smaller angle between R and dX 


• As for the direction of H, Biot and Savart observed this to be at right 
angles to both d\ and R and determined by the advance of a right-hand 
screw that rotates d\ onto R through the smaller of the two possible 
angles between them. From the definition of the cross product, we can 
therefore write the Biot-Savart law as 


<fH(r) 


7»r x R 

4tt7? 3 


(8.5) 


Note that in the example of the diagram above, dH at P is directed 
outward from the page. You can always quickly check the direction 
of magnetic field from the following “right-hand rule”. If you imagine 
gripping the conductor with your thumb in the direction of the current, 
your fingers will all point in the direction of the magnetic field. 


• The total magnetic flux intensity at P, arising from all of the current 
elements along a path (usually a closed path) is therefore 


H 


W = / 


l'{r')dX x R r l'(r')dX x e R 


4irR 3 


= f 


4irlfl 2 


( 8 . 6 ) 


where e R denotes the unit vector from each current element to the 
point, P , at which H is to be determined. 
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8.4 Magnetization and magnetic susceptibility 

• You should notice the following correspondences, that will make our 
later work with static magnetic fields simpler, using the experience 
that we have already gained in electrostatics: 

H «-► E (8.7) 

B D (8.8) 

p e (8-9) 

• The medium surrounding the current does not affect the magnetic field 
intensity, H, but it appears in expressions for the magnetic flux den¬ 
sity, B. Contrast this with the electrostatic case. The reason for this 
is that the medium consists of atoms whose electrons orbit and spin 
within their shells: these movements of the electrons can be regarded 
as currents in their own right, and such currents will produce a field of 
their own. This field is known as the magnetization , M, and it adds to 
the effect of the H-field that is due to the main current of interest. 

• In a similar development to what we saw with electric polarisation , we 
can now relate B and H as if there were only free space in the region 
of interest, but with the magnetization adding to the H-field: 

B = /z 0 (H + M) (8.10) 

We can then write 

B = pH = p 0 p r ll = n 0 (l + Xm)H (8.11) 

where Xm is the magnetic susceptibility of the material of the medium, 
and the relative permeability of the medium, // r , is merely 1 + y m . 

• We will not deal further here with the physical causes of these magnetic 
effects in materials. You should, however, try to read something further 
in this area, because of its relevance to such interesting topics as data 
storage, superconductivity, plasma containment, NMR, etc. 
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e.g.l Magnetic field on the axis of a circular 

current loop 


• Let the loop carry current 1 and have the 
z-axis as its axis. Its radius is a, and we 
consider the magnetic field at an arbitrary 
point P on the axis. Let P be distance 

2 from the centre of the loop. 

• In cylindrical co-ordinates, a current element 
on the loop has length dl = ad(f> and is 
always at right angles to the vector R that 
joins d\' to P. By Pythagoras’ Theorem, 

R has length y/a 2 + z 2 . Hence for each element 
of current, the magnitude flfdhe field intensity is 


H ' 





dH = 


ldd(p'R\ sin(90°) Iadcp' 


4'irR 3 


47 r (a 2 + z 2 ) 


( 8 , 12 ) 


• Integrating around the loop, all of the radial contributions to H cancel, 
so only the z-components are relevant. We have 

, TT . \ adH Ia 2 do 

dH : = dH sin 6 = , —. = = ---- (8.13) 

2 \ 9 


J[a 2 + z 2 ) 47r(a 2 + z 2 ) 


• To integrate for all values of <j >, we just multiply by 27r, yielding 

la 2 „ la 2 


H. = 


2 (a 2 + r 2 ■ 


or H = 


2 (a 2 +~ 2 )t 


(8.14) 


Note that in the centre of the circular loop (where z = 0), the H-field 
therefore has magnitude 7/2a, and is entirely in the z-direction. This 
is in accordance with our expectation from the “right-hand rule” 
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e.g.2 Magnetic field of a straight line of current 


• The geometry suggests the use of cylindrical co-ordinates. 


• We consider a point P, distance r from the centre of the straight line. 
With I flowing as shown below, H will come out of the page by the 
“right-hand rule”. There is no z-component of H since P is opposite the 
mid-point of the line, with the result that all z-components cancel. H 
is perpendicular to the page, and therefore has no r-component either. 
Hence, by equations 8.4 and 8.6. 


H* 



I sin 4 > dz 
AtvR 2 


(8.15) 


Then, since sin 4> = r/R and R = yfr 2 + z 2 , this gives 
Ir f a dz I 1 


H„ 


Jr r a 
47r J-a 


(r 2 + z 2 )i 27rr v ^y^j2 + T 


(8,16) 


from integration tables. 


• This is the expression for the magnetic field of a finite line of current. 
In the case of an infinite line of current, a —* oo, and then Hj, becomes 


n* 1 " 1 ~ 2 -r ° r H “ 


(8.17) 


T 



-a 




Chapter 9 

Ampere’s Circuital Law 

IV r ' 

9.1 The line integral of magnetic field 


• As a consequence of the work done by Ampere, it can be shown that 
the line integral of static magnetic field, taken about any closed path, 
must equal the current that the path encircles. Thus: 

j> H ■ d\ = I = ^ 3 ■ dS (9.1) 

• Notice that I is positive if the line integral is performed clockwise , 
looking in the direction of I. Notice also the similarity between this law 
and Gauss’s Law in electrostatics: in highly-symmetrical geometries 
this law. know as Ampere's circuital law , can be used to obtain the 
magnetic field very easily. 

• The proof of this law is not given here, as it must wait for a formal 
understanding of magnetic vector potential. 
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e.g.l The magnetic field due to a line current 


• Already calculated in the previous lecture, we now find this magnetic 
field by integrating around a circular path of radius r centred on the 
axis of the wire, as drawn above. Cylindrical co-ordinates is the natural 
system to use. The symmetry dictates that H will not vary with the 
angle d as we integrate around the path, and that H || d\ everywhere 
along the circular path. Hence, the integral becomes simply 

f H • d\ = H* J dl = 2?r rH* = I (9.2) 

from which Hs =-and Ba, — ——, just as we had before. 

2 nr 2irr 


e.g.2 Magnetic field between coaxial cylinders 

• Consider a coaxial cable in which current 1 
travels in the inner cylinder and returns in 
the outer. Exactly the same development as 
above, integrating along a circular path that 
lies between the cylinders, will give 

H 4, = -J— for a < r < b (9.3) 

2 nr 

• Note that, if r > 6, then the total enclosed current is 1 — I = 0, so that 
f H • d\ = 0, from which 



H = B = 0 for r > b (9-4) 

This explains why coaxial cable is used to transmit any signal from 
which components or systems that lie near the outer sheath of the 
transmission line need to be electromagnetically shielded. 
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e.g.3 The magnetic field of a solenoid 


• If current I flows in an infinitely-long coil (known as a solenoid ), which 
is wound with n turns per metre, then it is acceptable to regard the cur¬ 
rent as flowing in a sheet , entirely circumferentially around the solenoid. 
We ignore the tiny component of current flow in the z-direction that is 
due to the progress of the solenoid winding. The sheet’s current den¬ 
sity is nl Am -1 , and so the current that flows in a differential strip of 
length dz metres is nl dz amperes. 



• By the methods of the last lecture, we can write the field due to this 
differential element (which is just a ring of current) as 


dH z 


nla 2 dz 
2 (a 2 + z 2 )^ 


(9.5) 


and so, if the solenoid has length 26 and is centred at the origin, 


H z 


- L 


6 nla 2 dz 


-b 2 (a 2 + z 2 )2 
nla 2 f b 

(a 2 + z 2 )i 
2 


2 

nla 2 


r b 1 

J-b (a 2 + Z 2 )l 


dz 


2 a 2 


L[(a/6) 2 + l]-J 


(9.6) 

(9.7) 

(9.8) 


a T 

J r 

-A-'" - 
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• The integral in the line above is relatively easy to obtain, using the 
substitution z = a tan 9. Having found the expression for H z for a 
finite solenoid, we now let |6| —* oo to arrive at a simple expression 
(purely in the z-direction) for the for the magnetic field of an infinite 
solenoid at any point on its axis: 

H z = nl on the axis (9-9) 


• The diagram shows a sectional view of a 
part of an infinite solenoid, with the current 
directions indicated by dots and crosses. 
Now if we consider the 1-metre-long path 
shown on the left-hand-side of this diagram, 
we can use Ampere’s circuital law to give us 





jftd\ = nl (9.10) 

But we already know from the above that H z = nl on the solenoid’s 
axis, and that since H is purely z-directed there are no contributions 
to the integral from the radial sections of the path. Thus, all of the 
contribution to /H-dl comes from the axial section of the chosen path, 
and hence the field outside the solenoid is zero. 


• Furthermore, if we took the other 1-metre-long path shown in the dia¬ 
gram, then by Ampere’s circuital law we know that /H • dl — 0 (since 
the path encloses no current). But we know that /H • dl = nl for 
the axial section, so the contribution from the other horizontal section 
must be — nl. In addition, we note that the choice of where the upper 
horizontal section is drawn is arbitrary. In other words, the magnetic 
field is constant everywhere inside the solenoid , and of value nl Am -1 . 




Chapter 10 
Inductance 

10.1 Definition of inductance 


• In circuit analysis, we would describe any structure which stores energy 
in a magnetic field as an inductor. Recall from basic physics how any 
current has a magnetic flux ($) associated with it. Where two (or more) 
conductors are separated by free space, or by a material medium, the 
magnetic fluxes due to each conductor are said to be linked. Thus, for 
a coil of N turns for example, the total flux linkage is written 

A = N$ (10.1) 


• Since the flux, $, is usually proportional to the current, 7, flowing in 
the inductor, so too is the total flux linkage, A. Thus we may write 

X = LI (10.2) 

where the constant of proportionality, L in henries , is known as the 
inductance of whatever arrangement of conductors is being studied. 
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• Now the total flux linkage, A, can be found by simply adding up the 
magnetic flux density over the area in which the flux linkage exists: 

\ = J s B-dS (10.3) 

So, where the surface S has been specified, we have the important 
expression for inductance: 

L = \ / B -dS (10.4) 

1 Js 


• Take, for example, the coil of N turns shown here. Almost all of the 
flux due to the current passes through the shaded surface S. A little, 
however, passes through the wires of the coil itself. We often distinguish 
these as giving rise, respectively, to the external inductance and to the 
internal inductance of the coil. As suggested in the case of the coil of 
N turns, the internal inductance may be very much smaller than the 
external inductance, and is sometimes negligible. 


S 
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e.g.l The inductance of a parallel-plate 

transmission line 




Consider a pair of parallel plates, such 
as. for instance, in a microstrip transmission 
line. Let the plates have width w and have 
separation d. The separating medium has 
permeability p. Each strip is of thickness 
h , considered negligible, and a current / 
flows in the one strip and returns in the other. 





We make the important assumption that the 
strips are wide enough, compared to their distance 
apart, that the magnetic field lines between them 
are parallel. We write the magnetic field strength 
between the strips as H a , and note that we have 
made the assumption that d << w. We also assume 
that the magnetic field is virtually zero elsewhere 
- that is, we neglect any fringing field. 


. x 

C- 


tzZZZZZZZ22$ 


/' $ZZZZZZZZ2&' 

* -w-» 


4 


4 i 

u 

1 r 


• Then, by Ampere ’s circuital law, 

I = j H ■ dl = H 0 w 


(10.5) 


over a path of total length w between the conductors, and hence 



w 


( 10 . 6 ) 


• Now. noting that B 0 is parallel to dS everywhere over the the one- 
metre-long flux-gauging surface shown, we can find the inductance per 
metre of transmission line: 

L = 7 / B 0 • dS = ^ / dS = d x 1 ) = —— H/m (10.7) 
1 Js I Js wl w 
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e.g.2 Inductance of a coaxial transmission line 


• We know from previous examples that the magnetic field in a coaxial 
cable is entirely circumferential, and of value 


H* 


I 

2n r 


so 


Bj, = 


tL 

2irr 


( 10 . 8 ) 


where r lies between the radii of the two conductors in which the cur¬ 
rent I flows, and // is the permeability of the medium between the 
conductors. 


• Referring to this diagram of a coaxial cable, the total flux linkage per 
metre is to be found by integrating the flux density over the area S and 
noting that B is parallel to dS everywhere on the surface. 



(10.9) 

( 10 . 10 ) 


• Note that this expression gives the external inductance of the coaxial 
cable, because it only considers the flux in the region between the two 
conductors. We make the further assumption (particularly valid at high 
frequencies, as we shall see) that the internal inductance is negligible. 







Chapter 11 

The Curl, and Stokes’ Theorem 


11.1 Using and interpreting curl 

• Ampere’s circuital law gives rise to useful expressions, in integral form, 
for the magnetic field vector. In electrostatics, we saw the value of 
also having compact expressions for unknown quantities in differential 
form, and we introduced the vector operation of divergence to achieve 
this. In a somewhat similar way, the curl is used to obtain differential 
expressions for the magnetic field vectors, H or B. 

• It is assumed that you have worked with curl F (written V x F) before, 
particularly in MAM280W. If this is not the case, then you are advised 
to consult mathematics texts to gain a better understanding of the curl 
operation. Note, in particular, that V x F is a vector. 

• The basic definition of curl F, at a point and in the x-direction, is the 
limit of the line integral of the field F taken around the point, as the 
area enclosed by the path, and perpendicular to x, vanishes to zero. 
(Compare our earlier definition of divergence with this). Thus 

[ cur l F ], = (11.1) 

and this can also be stated in the y- and z- directions, to give 

curl F = V x F = [curl F] r + [curl F] y + [curl F] z (11.2) 
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Equivalently, although we will not prove this here, 


curl F = V x F = 


x 

a_ 

dx 


F x 


y 

a_ 

dy 

F, 



(11,3) 


• WARNING! Unlike divergence, the curl of a vector field does not have 
an easily-visualised interpretation. Curved field lines may possess zero 
curl, while some fields of straight, parallel vectors may have non-zero 
curl. Perhaps the best way to perceive the curl of a vector field at a 
point is to imagine that all of the field lines are velocity vectors and 
that a small, circular paddle sits at the point of interest. If the paddle 
spins, then the curl at that point is non-zero. 




As you might expect, there are equivalent expressions for V x F in both 
cylindrical and spherical co-ordinates, and these are given in formula 
sheets and (for V x H) below: 


VxH = 


?' sin 6 


£ { Htsm6)-^]+i[^^-£(rHt)]+$[£(rH e )-2&] (sph) 


VxH 


; I 
r 


QH Z 

d<p 



+ 4 >[ 


irafft 

dz 


8H z 

dr 


+Z 


1 1 dHr 
r dr r d<t> 


(cyl) 


(11.4) 


11.2 The curl of electrostatic field 


• Any vector, F, which can be written as the gradient of a scalar function, 
£(x,y,z), has zero curl. This is relatively easy to show, since, if 


_ *d( A <9£ „ d( 

F = v? = x & +y aJ + 2 & 

then, by our determinant definition of curl, 

( 


(11.5) 


VxF = x 


\ . (in a SL\ - (AL_ AL\ 

J ^ \dzdx dxdz) \dxdy dydx) 

( 11 . 6 ) 

o2 / o2 t 

But since, for a general scalar function /(a, b ), the equation 


\dydz dzdy 
or a g( 

above becomes simply 


V x F = 0 


(11,7) 
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• Now, electrostatic field is the (negative) gradient of a scalar function, 
known to us as the potential function, ${x,y,z). Thus, knowing that 
E = —V<f>, we have 

V x E = 0 (11.8) 

• You may recall from mathematics that zero curl is a characteristic of a 
conservative field. It is now easy to see why. Since 

<f E • d\ 

V x E = lim —- = 0 (11.9) 

as— ►() AS 

it follows that f E - dl = 0 , which is the condition for a conservative 
field. 


11.3 The curl of magnetostatic field 

• From the basic definition of curl, we see that the x-component of the 
curl of magnetic field intensity is 

[curl H] x = lim £5-^ (11.10) 

1 J as*—> o A S x ’ 

• But we saw in stating Ampere’s circuital law how § H • dl is related to 
the current density vector J = (J x ,J y ,J z ) by 

ji H ■ dl = J^3 ■ dS = 1 (11.11) 

Thus, 

i curlH i*=Je. i ^r- J * (1U2) 

by definition, and likewise [curlHjj, = J y and [curlH] z = J z 

• Therefore, by equation 11.2, we arrive at the compact and useful con¬ 
clusion that, for magnetostatic fields, 


VxH = J or V x B = gJ 


(11.13) 
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11.4 Stokes’ Theorem 

• Conversion betwen the integral and differential forms of relationships 
for electric fields was made easier by the divergence theorem. In a 
similar way, Stokes’ Theorem gives us a convenient “short cut” in our 
vector algebra when handling magnetic fields. 


• Imagine a surface, S, embedded in 
a field, F. We wish to divide 
the surface into infinitesimal vector 
areas, dS, and to calculate 
(V x F) • dS in each area. 


• Since, by the definition of curl 

(V x F) • dS —► <f F • dl as o?S —* 0 
JdS 



(11.14) 


we see that (V x F) • dS can be found in each little area simply by 
integrating F around its small perimeter. 

If we now add up all the contributions of (V x F) • dS for all of the 
small areas, then we find that at every boundary between small areas 
the line integrals cancel , leaving only the integral of F along the outer 
boundary of the surface. Hence 


J (V X F) ■ dS = j F ■ <9 


(11.15) 


• This is Stokes ’ Theorem. We can now apply it to magnetic fields and 
write 

j>n-d\ = j s {W xH)-dS (11.16) 

and, since V x H = J, we arrive back at a statement of Ampere’s 
circuital law: 

fH-dl = J3-dS = l 


(11.17) 
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e.g. Find the curl of the magnetic field both 
outside and inside a long, straight wire of 
radius a, that carries a current / 


• Outside the wire, the situation is very straightforward. The magnetic 
field is circumferential only and we know that H# = or equivalently, 
that H = We can find V x H in cylindrical co-ordinates by 

observing that H has no z- or ?- component and no z- or ^-dependency. 

dH z dH, 8H r dH z dH r n /iiioN 

(1L18) 


Also, since rH 0 = ^,we obtain 


d(rH<,) 


(11.19) 


Hence, outside the wire all components of the curl are zero, so 

V x H = 0 (outside) 


( 11 . 20 ) 


Inside the wire, there is a uniform current density J = ^ acting in 
the z-direction. Considering a circular loop of radius r < a within the 
wire, together with Ampere’s circuital law, reveals that 


H • d\ = 2'KrH r i > - L nr = 7 rr 2 |J| = —I 


( 11 . 21 ) 


Hence = 5 ^ 7 , which has ^-component only, and no z-dependency. 
Thus, in cylindrical co-ordinates, 

V x H = z d ±M _ i«*l = *1 * [All = (A) * ,11.22) 

r or r ocp J r or [27rorJ \7raV 
and so, as we should expect from our earlier equations, 


V xH 


z = J (inside) 


(11.23) 
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Magnetic Potential 


12.1 Defining vector magnetic potential 


• Just as we were able to calculate electrostatic fields knowing the elec¬ 
tric potential function $ and using E = — V$, so it is useful to define a 
magnetic potential , A, from which magnetic flux density can be calcu¬ 
lated. Magnetic potential is a vector , which arises because of a current 
flowing: each element of the vector A has the direction of the small 
element of current path that produces it, and we can evaluate A at a 
point by integrating the contribution from all the elements of the path. 


• We start with Ampere’s law 


H 


< r > = / 


/'(r')dl'xR 
47 tR 3 


or 


B 


< r ) = "/ 


/'(rQdl'xR 
47 tR 3 


( 12 . 1 ) 


which expresses the field intensity at a point r due to a current I flowing 
in a loop /, with each element of the loop at r' denoted dY. Recall that 
the vector R is just 


R = r - r' = x(x - x') + y (y - y') + i(z - z) (12.2) 


Then you can use the definition of the gradient operator to show that, 
in rectangular co-ordinates (try this on paper), 


dl'xR 

R 3 


= V 



xdl' 
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(12.3) 
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• Now, for any vector F and scalar function ip in rectangular co-ordinates, 
you can also show (and should do this for yourself) that 

Vx(^F) = ipV x F - F x Wip (12.4) 

Hence, replacing ip with ^ and F with dl', we have that 

Vx^ = -j-Vxdl'-dl'xV (4) (12.5) 

R R \RJ 

so, recalling that for vectors Vi and V 2 , Vi x V 2 = — V2 x Vi, we get 

v (s) x<n ' = Vx (f)-s Vx<i1 ' (12 - 6) 


• Back in Ampere’s law, this means that we can replace 


dl'xR 

R 3 


with 


V 



xdl' = Vx 



4vxdl' = Vx 
R 



(12.7) 


since V x dl' = 0 (because dl' is not a function of x, y or z and so has 
partial derivatives of zero with respect to each of these variables, so 
that its curl is zero too in rectangular co-ordinates). 

• Hence, Ampere’s law becomes: 


B(r) 



= V x A(r) 


( 12 . 8 ) 


if we define 




Z'(r') dl' 
R 


(12.9) 


We call A the vector potential of the magnetic field. In the case where 
we only know the current density J through the source of magnetic 
field, since I = J dS where dS _L J, and since dl' || J, dS and dl form 
a volume element dV , we can thus also write 



J(r') dV' 
R 


( 12 . 10 ) 
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12.2 Using B = V x A 

• Visualising magnetic vector potential is extremely difficult to do. In¬ 
stead, we just work with A as a useful intermediate function and re¬ 
member that B = V x A. This may make our work in calculating B 
actually easier. 

• Thus, for example, to calculate the magnetic field due to a single ele¬ 
ment of current /, flowing in the ^-direction, we could first obtain 

i/N f Lil'iv') dY A „ , 



We then observe that A is in the ^-direction but is dependent only on 
radial distance r. To obtain V x A, we will need to get A expressed 
in a suitable co-ordinate system, which would be spherical co-ordinates 
in which A r = A z cos6 and Ag = — A z sm6. Note that 6 is the angle 
that the vector from the source to the point of interest makes with the 
direction of the current in the source. Note also that A# = 0 since the 
current element has circular symmetry in the direction. 


• Using the definition of V x A in spherical co-ordinates, we therefore 
have ^ « t,: <" • f' 


<t> 

B = V x A = - 


d dA r 

dr rAc ' ] !M 


= (**•«) 


But this is precisely \dz\ |i?| sin 6, acting perpendicularly to both 

z and r. as you would obtain from applying Ampere’s law directly. This 
illustrates that B = V x A can be used as a possible alternative route 
for calculating magnetic field. 
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12.3 The divergence of magnetic flux density 


• We have already shown that, for a vector function F and scalar function 
£, where F = V£, V x F = 0. 

• Consider now Stokes’s theorem for a 
surface S bounded by a short closed 
line, as shown. Suppose the surface 
to be immersed in a field F. Then 


f F • d\ = x F) • dS 



(12.13) 


As we let the boundary line shrink to zero, the surface S becomes a 
closed surface, and the line integral shrinks to zero. Hence 


jf (V x F) • dS = 0 


(12.14) 


regardless of the nature of F or of the closed surface S. We can now 
apply the divergence theorem for any vector G 


i GdS =l 


V-G dV 


(12.15) 


to obtain, with G = V x F, 

j^(VxF)-dS = ^V-VxFdV = 0 (12.16) 

Hence, in this completely general case we have established that 

V • V x F = 0 (12.17) 

for any vector field F. 


• Now, since B = V x A, we can therefore write 

V • V xA = VB = 0 (12.18) 

This extremely important relationship tells us that there are no sources 
of magnetic flux, and therefore nothing which corresponds in magne¬ 
tostatics to the phenomenon of charge in electrostatics. Put another 
way, V • B = 0 is a statement that magnetic field lines will always form 
closed loops. 
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e.g. From the vector potential, find the magnetic 
field due to a uniform, circularly-symmetric 
axial current density 


• We are considering, in effect, the field inside a 
wire due to a current 1 flowing in the wire. Let 
the wire have radius a, and let the current density 
be J 0 = acting in the z-direction. We have 
already calculated H in such a wire, using 
Ampere’s circuital law, but here we will obtain 
the same result by means of calculating the vector 
potential A. 



• Now, it can be shown that the magnetic vector potential inside the 
wire obeys Poisson’s equation, V 2 A = —fiJ, and that in this case A is 
purely z-directed and of value 

A, = ^(x* + y*) (12.19) 

(You can check that V 2 A = —fx J in this case.) Then 

H = —B = — (V x A) (12.20) 


1 


x y 

JL JL 

dx By 

0 0 


z 

B_ 

dz 


Hjfk ( x 2 + y 2 ) 


( 12 . 21 ) 



( 12 . 22 ) 


-Jo 

2 


(xy - yx) 


(12.23) 
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• Converting this result to cylindrical co-ordinates, in which x = r cos <f>, 
y = r sin <j> and \ (f|x + ||y) = sin <^x + cos 4>y we obtain 


xy — t/x = r<j> 


and so 



Ir 

2na 2 


<t> 


exactly as we had before. 


(12.24) 

(12.25) 


12.4 Scalar magnetic potential 

• It is sometimes convenient to work with another quantity, known as 
the scalar magnetic potential , and written $ m , by analogy with scalar 
electric potential <h e . In current-free regions we can write (as we have 
seen) 

V x H = 0 (12.26) 

• We have also seen that any vector that can be expressed as the gradient 
of a scalar function has zero curl. The reverse is also true, and so we 
can write H as the gradient of a scalar function, and this is 4> m . Thus 

H = (-)V$ m (12.27) 


• Then, since V • B = 0, we can write 


V ■ B = V-//H = V- (-^V$ m ) 

(12.28) 

and so 


V 2 $ m = 0 

(12.29) 


from which we can also see that $ m , like $ e , obeys Laplace’s equation. 
There are methods available (as we saw in electrostatics) for solving 
Laplce’s equation, so yet another way of calculating magnetic field may 
be to obtain in a given situation, and then to find H = — V$ m . 




Chapter 13 

Energy in Magnetic Fields 


13.1 Energy storage in static magnetic fields 

• You may recall that we derived this expression for the energy con¬ 
tained in a small volume of a region in which there is an electrostatic 
field: cIUe = |D • E dV. One can reason in a directly analogous way 
with static magnetic fields, provided that B and H are proportionately 
related—that is, provided that /x is constant for the material medium, 
in which case the medium is described as being linear. We may then 
write B = /xH where /x is a constant and, with this proviso, we have 
(after integration and if B || H): 

U H =\jB-HiV=tj v HUV ( 13 . 1 ) 

• Many magnetic materials are not linear—in other words, 
the value of /.i depends upon the strength of the magnetic ® 
field applied. The result is that, if you plot |B| against ^ 

|H| as the field strength varies in this material, you will t 
obtain a curve like the one shown. As we shall see, this 
complicates the picture of the energy stored in the field. 

In the non-linear case, the extra energy stored in a 
volume V for a small increase in B is 

dU H = J v B.-dBdV 
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( 13 . 2 ) 



j w 
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• If H || B, which is a feature of isotropic materials, we can then write 

dUn = H dB J dV (13.3) 

- /\r n . 

which means that the shaded area in the graph above represents the 
extra energy stored per unit volume for a small increment in the B-field. 


13.2 Energy loss in hysteretic materials 

• We have therefore seen that a non-linear, isotropic material will respond 
to an increase in field strength from zero by travelling from a to 6 in 
the B-H graph shown below. 



• By the time the material is brought to point 6, small magnetic domains 
within the material may have been aligned by the applied field. There 
is a tendency for this alignment to be semi-permanent, with the result 
that the material is somewhat magnetised , i.e. it will act as a permanent 
magnet. 

• If the applied field is now decreased, the material does not retrace the 
path from b to a, but will travel the path cde as the applied field is 
brought first to zero and is then made negative. Note that at point 
c. where the applied H-field is zero, there nevertheless remains a 13 - 
field due to the magnetisation of the material. This is known as the 
remanence or retentivity of the material. It then requires a further 
applied field in the negative direction to bring the B-field to zero at 
point d. This value of H is known as the coercive force. 
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• Similar arguments apply by symmetry as the applied H-field is again 
increased, and it is found that the material will trace a path from e 
through / and g back to b. The tendency of the material to “remember” 
its most recent state of magnetisation and not to retrace its B-H path 
until it has traversed the full loop bcdefg is known as hysteresis , and 
the loop is called a hysteresis loop. 

• Applying our law for energy storage as the B-field varies, we see, for 
example, that as the material moves from g to b there is a point (shown 
on the figure) where the increase in B-field leads to the storage of energy 
H\ dB per unit volume. Then, as the material moves from b to c, at 
the same B-field strength there is a release of stored energy equal to 
H '2 dB. The difference between the energy added to the stored energy 
and the energy taken from storage is (Hi — H 2 )dB, which is shown 
in cross-hatched shading on the diagram. This energy is effectively 
wasted, being dissipated as heat during the hysteresis process. 

• It should therefore be clear that ii, one transit of the hysteresis loop 
bcdt fgb , the energy dissipated per unit volume of the material is equal to 
the area enclosed by the hysteresis loop. Materials with narrow hystere¬ 
sis loops are termed “soft”, and are of value in electromagnets, trans¬ 
formers and motors because their energy wastage is relatively small. 
“Hard” magnetic materials have broad hysteresis loops and are used 
for making permanent magnets. 


13.3 Magnetic field energy and inductance 

• In the analysis of linear circuits, we are familiar with the expression 
\LP for the energy stored in the magnetic field of a circuit element. 
Since we have already seen how to calculate the inductance of certain 
symmetrical structures directly from field considerations, we now have 
the opportunity to check that our two definitions for stored energy are 
the same. 
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• The examples which follow will demonstrate and explore the useful 
relationship: 

[ —H 2 dV = -LI 2 (13.4) 

Jv 2 2 

• Note carefully that the inductance L refers to the external inductance 
of the structure if we consider only the field H that lies outside the 
conductors, and refers to the internal inductance if we consider only 
the H-field within the conductors. Because the external field generally 
occupies a much larger volume than the internal field, the external in¬ 
ductance usually far outweighs the internal inductance. The equation 
above applies equally to both situations, however, and, because of lin¬ 
earity, will also apply for the total inductance , if the whole volume both 
within and outside the conductors is considered. 

e.g.l Show that the relation fy dV = |L/ 2 
holds for the region between the two 
conductors in a coaxial cable 

• We already know from Ampere’s circuital law that the external induc¬ 
tance is L = 57 In H/m and that the field between the conductors is 
H$ = 577 - For a 1-metre-long cylinder V = nr 2 xl, and so dV = 2irr dr. 

• Hence, if we consider a 1 -metre length of coaxial cable which has an 
inner conductor of radius a and outer conductor of inner radius 6 : 

Utfdv = 

Jv 2 


✓T, 

\c 


\aV_ r 1 

2 Jv (2nr) 2 


iv = tpf^L. dr 
2 J a (2irr) 2 


1X1 Ja 7 dr = Y X 27r ^ ln ^ “ ln 


2 x 27r 


p »i ( b ' 

J2i ]n {a 


\ LP 


(13,5) 
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e.g.2 Find the inductance per unit length of an 
infinite solenoid of radius a which has n turns 
per metre and is filled with a material 
of permeability //, 


• We recall that the field is H z = nl everywhere inside an infinite 
solenoid. The volume per unit length is V = 7r a? x 1. Hence for 
any one-metre length we can write the stored energy as 

I ^-H 2 dV = tL n 2 I 2 V = ^ n 2 I 2 Tra 2 (13.6) 

•J Li Lt Ld 

• But the energy stored is also \LI 2 , so 

-LI 2 = —n 2 I 2 Ta 2 (13.7) 

2 2 

from which 

L = fm 2 fta 2 = pin 2 A [H/m] (13.8) 


where A is the cross-sectional area of the solenoid. 


• Note that the inductance of a coil is related to the square of the number 
of turns per metre. 
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Chapter 14 

Time-Varying Electric &; 
Magnetic Fields 

14.1 Summary of important relationships 


• Up to this point, our analysis of fields has included no dynamic effects. 
In all of our examples, the sources of electric or magnetic field have 
fallen into one of the following categories: 

1. They have been stationary in space and invariant in time. 

2. They have varied in time very slowly (or at very low frequency), 
so that over a given period of interest they have been effectively 
static. 

3. They have moved in space, but, as with a steady current in a 
conductor, this movement has led to no change in the overall field. 
The field in this situation may be described as quasistatic. 


• Under static or quasistatic conditions, we have identified or derived a 
number of important definitions and relationships between field quan¬ 
tities. These are summarised on the next page: 
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Major Definitions 

1. Electric field: E = ^ N/C 

2. Conduction current: J = crE A/m 2 
(in a conductor) 

3. Convection current: J = p\ p A/m 2 
(in space or a medium) 

4. Permittivity: e = e r e 0 = §■ F/m 

(We assume e to be a constant, because we assume the material 
medium to be homogeneous, isotropic, linear and time-invariant.) 

5. Permeability: p = p T p 0 = g H/m 

(Although this is often not the case in practice, here we also fre¬ 
quently assume that the material medium is homogeneous, isotropic 
linear and time-invariant.) 


Vector relationships for static fields 

1. V D = p (encapsulates Gauss’s law and our understanding 

of charge and flux); 

2. V • B = 0 (expresses the absence of magnetic field sources 

that would be analogous to electric charge); 

3. V x E = 0 (describes the conservative nature of electrostatic 

field); 

4. V x H = J (the differential form of Ampere’s circuital law, 

relating magnetic field and current density). 
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14.2 Time variation in fields 

• It is an observed fact, first noted by Faraday, that a time-varying mag¬ 
netic field will generate an electric field. There is also a complementary 
effect by which a time-varying electric field gives rise to a magnetic field, 
which was first quantified by Maxwell. 

• It is easy to show that changing magnetic field produces an electric field. 
You can move a conductor in the vicinity of a magnet and measure the 
voltage induced in it using an oscilloscope. Alternatively, you can put 
the conductor near an electromagnet and get the same effect by varying 
the strength of the electromagnet. It is less easy to show how magnetic 
field results from a changing electric field, but we will shortly see that 
this is indeed the case, and we will be able to quantify the effect. 

• Our static equations make no mention of the interdependence of chang¬ 
ing electric and magnetic fields. We will now develop some ways of 
expressing this interdependency, with the aim of refining our vector 
relationships so that they take the effects of time-varying fields into 
account. 


14.3 Faraday’s law of electromagnetic induction 


• Faraday found that the voltage induced in a conductor that is linked 
by a time-varying magnetic field is 


V = 


d<H m 

dt 


(14.1) 


where \k m is the magnetic flux linking the conductor at any instant. 
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• In the case of a coiled conductor with n turns in the coil, the inductive 
effect is additive , so we write 

(14.2) 

Note that the magnetic flux that links the coil may change if 

— either a nearby source of magnetic flux (e.g. an electromagnet) 
varies in strength or position, 

— or the coiled conductor is moved in the magnetic field. 

• Note also that it is not necessary to have a conductor for a change 
in magnetic field to produce a voltage. The same voltage would be 
induced between points in an insulator, or even in free space, by the 
same change in magnetic field. 


V = n 


d^rr 

dt 


• When there is a contribution to voltages in a region due to changing 
magnetic flux, it is no longer true that the line integral of electric field is 
independent of the path (i.e. the electric field is no longer conservative). 
Instead, we define, for any closed path, the electromotive force which 
is the extra induced voltage: 

/f)\u 

E ■ d\ = (14.3) 




• Then, since the flux linking a loop, M/ m , is simply the integral of the 
normal components of flux density through any surface bounded by the 
loop, we can write 



/ 


E ■ d\ = 


dV r 


.-It 

dt Js 


B • dS 


dt 


(14.4) 
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• Applying Stokes’s theorem to the left-hand term above, we get 

/ s (VxE). rfS = -/ s f -<& (14.5) 

and, for this to be true for any surface 5, we can set equal the quantities 
under the integrals, and write 

V x E = Uc'A (14.6) 

• This supports our view of electrostatic fields as conservative because it 
confirms that V x E = 0 (if ^ = 0). However, in this more general 
expression, we see clearly that if nearby magnetic field is time-varying, 
then the line integral of electric field around a closed path is non-zero — 
this is because some energy is stored in (or delivered by) the magnetic 
field as it changes in strength. Thus, net work would be done on or by 
a charge moving around such a closed path. 



CHAPTER 14. TIME-VARYING ELECTRIC k MAGNETIC FIELDS 89 


e.g. Find the greatest rate at which the magnetic 
field of this large electromagnet can be varied 
from 0 to 10 T without breakdown of 

the surrounding air occurring 


B 


• The electromagnet has circular poles of radius 20 cm, and we assume 
that the magnetic field between the poles is uniform. The field else¬ 
where is taken to be zero. All magnetic flux thus passes at right angles 
through the circular ends of each pole. 



• For any loop of radius r between the poles as shown, Faraday’s law 

d 


f Ed '=~mJs Bds 


(14.7) 


gives 


2nrE ( h = —nr 2 


dB z 


nr 2 B„ 


at - r < 14 ' 8 > 

where B max = 10 T and r is the time taken to raise the field to that 
value from zero. Note that, since oc r in this expression, the greatest 
£’-field will occur at the greatest possible value of r between the poles, 
i.e. at r = 20 cm. 


• Since air breaks down at about 3000 kV/m, we can now write 

nr B max ^Brnax 0.2 X 10 1 , 

r “ 2nr\Et\ = 2|£*| = 2 x 3 x 10 6 = 3 ^ ^ 4 ' 9 ^ 

This is the shortest time in which the magnetic field can be varied from 
0 to 10 T without creating an electric arc in the region of the poles at 
r = 20 cm. 



















Chapter 15 


Continuity and Displacement 
Current 


15.1 The equation of continuity 



• A fundamental principle, that is never seriously challenged, is that of 
continuity. Since current is comprised of moving charge, if a current 
flows out of a region (through the region’s bounding surface), then the 
positive charge inside the region must decrease. Thus 

f dQ d r ' ‘ ■ 

£3 ■ dS = I = = —— j pdV (15.1) 

11*1 C 


Dividing through by the region’s volume, AF, and letting —»• 0, 

we get 


lim 

AV'-O 


f3-dS 

AV 


lim 
AV—*0 


( dfpdV\ 
\ dt AV J 


(15.2) 


• This is recognised, through the definition of divergence, to be the same 
as the equation of continuity. 

= (15.3) 

in which p is the net charge density (not just the density of mobile 
charge) in the small region of interest. 
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• Note that in conductors dp/dt is non-zero only for very brief transient 
periods. This can be seen by imagining that a solid conductor has a 
net charge density of p = p 0 at some time t = 0 in some region of itself. 
Then, since J = crE = crD/e, the equation of continuity gives 

V-J = V--D = —^ (15.4) 

e at 

Since a and t are constants, 
this leads to 

and hence 


and recalling the relationship V • D = p, 
ap dp 


D) = 


dp a 

m + l p = 0 


dt 


(15.5) 

(15.6) 


• This is a familiar first order linear differential equation in standard 
form, for which the solution is 

p = Po <rW* (15.7) 

and so we see that the net charge density, p, will decay from its initial 
value of p 0 exponentially , and with a time constant or relaxation time 
of e/cr seconds. This is the time after which only 36.8% of the initial 
net charge density remains: in conductors it is very brief. For example, 
the relaxation time for silver is e/c r ~ 1.44 x 10 -19 s. Clearly, after a 
very short time (e.g. 10 -17 s), both p and dp/dt will be virtually zero. 

• Thus, for good conductors only, we may write 


V-J«0 (15.8) 

You may recognise this as the field equivalent of Kirchhoff’s Current 
Law, which states that the net current leaving a node where several 
perfect conductors meet is equal to zero. 



CHAPTER 15. CONTINUITY AND DISPLACEMENT CURRENT 92 

15.2 The definition of displacement current 


• We know that, for a static magnetic field H in an arbitrary medium, 
we can write 

J c = V x H ' fv > (15.9) 

where J c is the current density due to the actual motion of charges in 
the medium. We call J c the conduction current density. Now, if the 
magnetic field were time-varying, then there would also be time-varying 
net charge density in the medium because of the induced currents that 
would participate in the electromagnetic variations (i.e. dp/dt ^ 0) 


• However, taking the divergence of equation 15.9, we see that 

V-J c = V-VxH = 0#-^ (15.10) 

Hence, the equation V x H = J c cannot apply to time-varying mag¬ 
netic fields, because it leads to a contradiction with the equation of 
continuity. 


• Maxwell postulated that time-varying magnetic fields cause the flow 
of another type of current, called displacement current , whose density, 
the displacement current density , he defined as 

(15.11) 

Maxwell imagined that this extra (unmeasurable) current density forms 
part of the response to a time-varying magnetic field, and so he amended 
the purely static relationship V x H = J c to read 

VxH = J c + J d (15.12) 



thereby generalising the equation to the dynamic case. 
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• Now. with 

J c = VxH-Jj (15.13) 

as Maxwell proposed, we can take the divergence of the equation and 
obtain 

V J C = V-(V x H)—V-Jrf = O-V-f = -j^(V-D) = (15.14) 

and so the new equation does not contradict the equation of continuity. 


• We can therefore replace the fourth static equation from the last lecture 
(V x H = J) with the improved relationship: 

VxH = J+® (15.15) 

bearing in mind that we revert to the static case whenever the electric 
field is not changing because we then have ^ = 0. Wherever we 
know that electric flux lines vary in a region over time, however, the 
displacement current term gains importance, and the new relationship 
should be used. 
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15.3 Maxwell’s equations 

• We are now in a position to re-list the four vector relationships that 
we established for static fields and stated in the last lecture. The third 
and fourth equations are modified to take time-variation of fields into 
account. These four equations, taken together, are always referred to as 
Maxwell’s Equations , and they summarise what is known about electro¬ 
magnetic fields. We will spend the rest of this course using Maxwell’s 
equations to gain a clearer understanding of certain important electro¬ 
magnetic phenomena. The equations are: 

f . r A r ?.f 

1. V • D = p 

2. V • B = 0 

3. V x E = —f 

4. VxH = J + f 



• Note that each of these equations may also be expressed in integral 
form, either by applying the divergence theorem or else by using Stokes’s 
theorem. You ought to practise performing all of these conversions in 
either direction. The integral forms of Maxwell’s equations are: 


1. § s D-dS = $ vP dV 


(Gauss’s Law) 


2. f s B-dS = 0 


(No Magnetic Charges) 


3. fE-dl = -|/ 5 B-dS 


(Faraday’s Law) 


4. /H-dl = / s J-dS + |/ 5 D-dS 


(Ampere’s Law) 


\7 K 


"J J S 
5 








DurW'^: 


y.T do 1 
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15.4 Visualising displacement current 

• The fourth equation above is a refined version of Ampere’s circuital 
law, applicable to the time-varying case. We will take a further look at 
it to try to get a better visualisation of the displacement current. 



• When an alternating source, v(t) = V^,sin(u;t), is applied to a capacitor 
as shown here, a time-varying conduction current, i c , is known to flow, 
where 

i c = = ljCV 0 cosM) (15.16) 

at 

and, since C = tA/d, this current has density of magnitude 

J c = ^ = u ~jV 0 cos(wf) (15.17) 

• Now, Ampere’s circuital law states that the integral of magnetic field 
around a closed path is equal to the current through any surface that 
is bounded by that path. Equivalently, it is equal to the normal com¬ 
ponent of current density, integrated over any surface bounded by the 
path. This should apply at any instant or over periods of time, and it 
must be completely independent of the surface chosen. 

• Clearly, if we choose to integrate only conduction current density over 
surfaces 5i and S 2 in the figure above, then we will get a non-zero result 
for Si in general, but zero for S 2 , since no conduction current can be 
flowing in the dielectric of the capacitor. 
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• We note, however, that the electric field in the capacitor, which has 
magnitude E = v(t)/d, is time-varying, so that 


dD 

dt 



(15.18) 


In fact, there is a displacement current within the capacitor whose 
density has magnitude 


_dD__ d_E_ 
Jd ~ dt ~ e dt 


e dv(t ) 
d dt 


LO-V 0 COs(ujt) (= J c ) 
d 


(15.19) 


Assuming that no fringing occurs outside the capacitor plates, the dis¬ 
placement current id may then be found from 


id = A|Jd| = u> K> cos (cot) = uCV 0 cos(u>t) (15.20) 


• But this is precisely the expression found for i c above! Therefore, con¬ 
tinuity of overall current is preserved at the capacitor’s plates. In the 
absence of fringing (so that D varies only inside the capacitor), and 
since J c = 0 within the capacitor’s dielectric, we see that ONLY i c 
contributes to the integral over Si and ONLY id to the integral over 
S 2 . This means that our requirement of consistency from Ampere’s 
Law is satisfied by the new version stated in the fourth equation above, 
because 


j)H- d\ 



d f 

dS + 0 = 0 + — / D dS 

dt Js 2 


(15.21) 


and, in this example, 



(15.22) 


• We come to an understanding, therefore, that displacement current is 
vital to the consistency of Maxwell’s equations, even though it is not a 
directly measurable quantity. 






Chapter 16 

Sinusoidally-Varying EM Fields 

16.1 The use of phasors 


• Time-variation of electromagnetic fields is very often either sinusoidal 
or else periodic, and in the latter case it is subject to Fourier decom¬ 
position into sinusoidal components. 

• You may recall from your circuit analysis course that the mathematical 
handling of sinusoidal functions is made easier by carrying the ampli¬ 
tude and phase information about the sinusoid in the form of a complex 
number known as a phasor. 

• Thus, the sinusoidal function v(t) = V m cos(u it -f <f>), which has ampli¬ 
tude V m and phase 4> with respect to some agreed reference, is repre¬ 
sented by the phasor 

V = V m L<f> - V m cos (<t>) +jV m sin(<t>) = V m e 34> (16.1) 
Note that the frequency, u>, is assumed known. 
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• The time-dependent function v(t), may then be recovered from the 
phasor representation by multiplying the phasor by e 3UJt and then taking 
the real part of the product. Thus 

v(t) = 


= ft(y m e j(u4+ ^) 

= 3£(Kn cos (cot + (f>) + jVm sin(u>f -f- cj))) 

= I'm cos (cot A <t>) (16.2) 

exactly as we defined it. 


• One motivation for working with complex exponential phasors rather 
than with sinusoids is that it makes the differentiation of field quantities 
with respect to time very much easier to handle compactly. The time- 
derivative of a phasor is simply the phasor multiplied by ju. In this 
way, 

dv(t) 

—— = -uV m sin (u>t A d>) (16.3) 


can be represented by ju\ because 


dv(t) 

dt 


= ${juVe ju,t ) 




= 3?( juV m cos(t ot A 4>) + j 2 u>V m sin(u4 + </>)) 


= -ioV m sin (u>t A (t>) (16-4) 

after the phasor representation is returned to the time-domain form. 
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16.2 Maxwell’s equations in phasor form 


• We can now re-write Maxwell’s equations for the case of sinusoidal 
time-variation. Note carefully that although we use the same symbols 
in these equations, they actually all denote complex numbers which 
carry the magnitude and phase information about each field quantity at 
a particular instant in time. To “unlock” the time-dependent function 
from each of these phasors, we would multiply by e ^ ujt and take the real 
part, as described above. 

• Maxwell’s equations then become, in phasor form: 

1. V D = p(x,y,z,t) 

2. V • B = 0 

3. V x E = —juj B 

4. V x H = J + jcot) 


• Note the compactness of this form of expression. Equation 3, for ex¬ 
ample, can be decomposed into 


V x E = 


x y z 

JL JL JL 

dx dy dz 

E x Ey E z 


(16.5) 


\ oy oz ) \ ox oz ) \ ox ay J 

= -juB x x - juBy y - juB z z (16.6) 


so that, for example, 


dE z 

dy 


dt 


dz 


- = -juB x 


(16.7) 


in which the quantities shown are still all phasors whose time-dependent 
equivalents can be unlocked in the usual way. 
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16.3 Plane waves in free space 

• Fortunately, it is often possible to make massive simplifications to these 
sets of equations and still to be able to study some very realistic situ¬ 
ations with them. 


• Take, for instance, the behaviour of an electromagnetic wave which 
varies along only one axis (the z-axis, say), and which exists in a region 
in which there is no free charge (p = 0) and also no movement of charge 
(J = 0). We call this a plane wave in free space. 


• It turns out that electromagnetic waves travelling in a vacuum (e.g. 
through outer space) are a good example of plane waves in free space. 
Since the interplay of electric and magnetic fields in such waves is gov¬ 
erned by Maxwell’s equations, we are motivated to apply Maxwell’s 
equations with p = 0, J = 0, = 0 and = 0 to see what can 

be said about light, microwaves, radio waves etc. when they travel 
through a vacuum (or near vacuum). 


• For plane waves in free space, Maxwell’s equations are given below. 
Their phasor equivalents, which we will use later, are given in brackets. 


1. V D = 0 

2. V B = 0 

3. VxE=—f = -//f 

4. V x H = f = cf 


(V • D = 0) 

(V-B = 0) 

(V x E = —ju>p, H) 
(V x H = jue E) 


• Now, remembering that d/dx = 0 and d/dy = 0 for a plane wave, the 
decomposition of the equations is a bit easier. Equation 3 now gives 


\ 


x y z 
0 0 & 
E x E y E z 




-\ n 

a ^ 





V x E = 


(16.8) 
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so that we have the three relations 


dE v 


-n- 


dH x 


d E x 


dz r dt ’ dz 

Similarly, equation 4 now yields 


= -fi- 


a f- o.-4Sr u«> 


V x H = 


x y z 

0 0 h 

H x Hy H z 


giving us three more equations: 




dHy dE x 


di 


= e- 


dt ’ 


dH x 

dz 


= e- 


dE± 

dt ’ 


= e 


dE 

dt 


(16.10) 


0 = e 


dE z 

dt 


(16.11) 


• We can now combine the second equation from V x E with the first 
from V x H, for example. Differentiating the second equation from 
V x E with respect to 2 gives 


d 2 E x d 2 H y 

dz 2 ^ dzdt 


(16.12) 


Likewise, differentiating the first equation from V x H with respect to 
time yields 


d 2 H y _ d 2 E x 
dzdt dt 2 


(16.13) 


We combine these results by eliminating H y to arrive at 

d 2 E x d 2 E x 
dz 2 “ //e dt 2 


(16.14) 


• This is a second-order partial differential equation, known as the one- 
dimensional wave equation , which we must solve if we wish to under¬ 
stand how E x will perform as part of the planar free-space wave. Note 
that we could just as easily have combined the above equations to ob¬ 
tain wave equations for E y , H x or H y . 
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16.4 Solving the 1-D wave equation 


• Let us suppose that some function F(z, t) satisfies the one-dimensional 
wave equation, so that we may write 


d 2 F _ 1 d 2 F 


(16.15) 


where v is a constant whose meaning and value are to be determined 
shortly. 


• Mathematical solutions to this equation take the form 

F(z,t) = f 1 (t- Z -)+f 2 (t + ^) (16.16) 

where ,/i is some function of (t — z/v) and f 2 is a function of (t + 
z/v), so all the brackets here denote arguments of the functions (not 
multiplication!). It is easy to see that this is a form that a solution could 
take. Simply re-substitute it into both sides of the wave equation (you 
should do this for yourself as a check). 


• Now think about j\(t — z/v), a function that varies in both space ( 2 ) 
and time (t). Imagine the value of this function at z = 0 and t = 0. A 
moment later, t has moved on, so (t — z/v) will have changed, and so 
the value of the function at z = 0 will now be different. If we wanted to 
“stay with” the same value that we had at z = 0 and t = 0, evidently 
we should have to keep (t — z/v) constant. But if 


then 



v 


(a constant) 


(16.17) 


z = vt — K v 


(16.18) 


and we see that we should have to move along the z-axis with velocity 
dz/dt = v in order to “stay with the same point” on the /1 function. 
This is tantamount to saying that f\(t — z/v) moves along the z-axis 
with velocity v. 


• By a very similar argument, f 2 (t+z/v) can be shown to move backwards 
along the z-axis, with velocity —v. Check that you can visualise this. 
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16.5 Velocity of electromagnetic propagation 


• We can now apply all of this theory to plane electromagnetic waves in 
free space, travelling in the z-direction. We know that, in this case, E x 
satisfies the wave equation 


d 2 E x _ d 2 E x ^ i d 2 E x 
dz 2 ~ ^ dD ~ V 2 dt 2 


(16.19) 


• Most commonly, electromagnetic waves have sinusoidal time-variation, 
so that the solution to the wave equation, which is of the form 

Ex = h (t ~ ") + h (i + -) (16.20) 

has sinusoidal components j\ and /Y 

• This suggests the use of phasors to solve the wave equation, which can 
be written for this purpose as the one-dimensional Helmholtz equation 

= -LU 2 pitE x (16.21) 

simply by replacing time differentiation with multiplication by ju>, and 
then recognising that E x is now a function of one variable only (z), so 
that the partial derivative can be replaced by a total derivative. 

• Because E x is a phasor in complex exponential form, the solution com¬ 
ponents fi and /2 should also be written as complex exponentials, and 
are also functions of z only. Hence 

E x = Cl e~ jkz + c 2 e ikz . (16.22) 

where k = w^f\xi. You should verify, by re-substitution, that this 
solution does indeed satisfy the Helmholtz equation. 

• We can now recover the time-varying solution for E x (t) from the phasor 
E x by simply writing 
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E x (z,t) = $(E x e jMt ) 

= n{ Cl e~ 3kz e ]Wt + c 2 e jkz e jwt ) 

= S(c 1 e J ' M_fa) + c 2 e j{ujt+kz) ) 

= Ci cos(u;f — kz) + C 2 cos (u>t + kz) (16.23) 


which is, exactly as expected, a pair of sinusoids, each of frequency u, 
and travelling in the +z and — z directions respectively. 

• The velocity of these waves is found, by comparison with the wave 
equation, to be v — 1 / y/ye, and in a vacuum this takes the value 

v = —L= (16.24) 

y/i 

_l_ 

V47T x 10- 7 x 8.854 x 10~ 12 

w 2.9979 x 10 8 [ms -1 ] 

We therefore have the very satisfactory result that light (or any other 
electromagnetic wave) travels at the velocity expected from field theory. 
We can also predict the velocity of electromagnetic propagation in any 
other medium, if we know the values of e = e r e 0 and y = y r y 0 in that 
medium. 


• Note, finally, that k, which is known as the angular wave number of 
the travelling waveform, signifies the number of complete waves in 2ir 
metres along the z-axis and can be written 

k = ujy/ye = — (16.25) 

v 

and that therefore the distance between similar points in a periodic 
wave, which is given by 27r//c, is 


2tt 2ttv v 

= T = ~ = 7 


(16.26) 


where / = u;/27r is known as the frequency and A is the wavelength of 
the periodic waveform. 









Chapter 17 


The Skin-Depth Phenomenon 


17.1 Maxwell’s equations in good conductors 

• The previous chapter dealt with Maxwell’s equations in free space—an 
excellent dielectric. We now study the behaviour of electromagnetic 
fields in the vicinity of a good conductor. The conductivity of such 
a medium is given at frequencies up to many gigahertz, by the real 
number a, (although it may be noted that a is complex at optical 
frequencies). We can therefore state the familiar field-theory version of 
Ohm’s Law in phasor form: 

J = <rE (17.1) 

• Now since V x H = J + ju>t E (the 4th of Maxwell’s equations), a 
convenient factorisation becomes possible, and so we can write 

V x H = (cr + jut) E (17.2) 

It is important to note that in good conductors cr toe in the bracketed 
term above, while in good dielectrics J is insignificant and we have been 
able to assume that ut 

• We have a picture of how conduction takes ^ ^ ^ ^ 

place in a conductor through the drifting 
charge carriers within a lattice of atoms 
when a field is applied across the lattice. 
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Carriers are promoted from a valence to a conduction role in the ma¬ 
terial by the addition of energy in the form of heat or light. The total 
number of carriers (electrons) in a given volume V equals the total 
number of positively-charged nucleons, however. As long as V is large 
with respect to individual atoms yet small in comparison with major 
circuit elements, the net charge within V will therefore be zero. This 
picture is in good agreement with our field theory, as can be seen by 
taking the divergence of both sides of the equation above: 

V-VxH = (<T + jwe)V-E (17.3) 

Since V • V x H = 0 by definition, we arrive at V • E = 0 and so 

V.D = p = 0 (17.4) 

confirming that the net charge, p, will be zero when Ohm’s Law applies. 

17.2 E, H and J in good conductors 

• It is quite straightforward to find differential equations which describe 
the behaviour of field quantities E, H and J in the vicinity of a large, 
good conductor. The general pattern for deriving these relationships is 
to work with Maxwell’s equations, Ohm’s Law and the vector identity 
V xV xF = V(V • F) — V 2 F that we met earlier in the course. You 
ought to gain confidence in all three of the derivations below. 

I: In the case of the electric field intensity, E, we can begin with the third 
Maxwell equation, V x E = —jw/xH, and then take the curl of both 
sides to get, 

V x V x E = V(V • E) — V 2 E = —yw/xV x H (17-5) 

Then, since Ohm’s Law leads to V • E = 0 (as we saw above) we obtain 

V 2 E = jw/x(cr + jwe)E (17-6) 

and you will notice that we have also replaced VxH with the Ohm’s 
Law version that we derived above. Now in a conductor a we, and 
so we can effectively neglect displacement current altogether and write 

V 2 E = jfw/xcrE (17.7) 






6 ( 4 ^ 6 ) 
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II: Similarly, for the magnetic field intensity, H, we might begin with 

V x H = (a + ju>e)E ~ <tE (17.8) 

Taking the curl of both sides now gives 

V x V x H = V(V- H) - V 2 H = <rV x E (17.9) 

But V x E = —j to/uH by the third Maxwell phasor equation, and 

V • H = 0 by the second, and so 

V 2 H = +jui fiall (17.10) 

III: Finally, for the conduction current density, J, a possible starting point 
is to take the curl of both sides of Ohm’s Law and then to apply the 
third Maxwell equation, V x E = — jufitl, to obtain: 

V x J = <tV x E = — ju/iail (17.11) 

Next, we can take the curl of both sides again, and this time replace 

V x H to get 

V x V x J = —jtofiaV x H = — ju>fia(a + ^u;e)E (17.12) 

Now, remembering that a ut for a conductor, we can neglect the 
displacement current term juie E and replace the other term, crE, with 
J (by Ohm’s Law) to get 

V x V x J = V(V • J) - V 2 J = -jujfiaJ (17.13) 

Then, finally, we must appeal to the concept of continuity , in which we 
saw that V • J = — |y w 0 for a good conductor in the steady state. 
This then leads us to 

V 2 J = jujfiai (17.14) 

• The differential equations governing E, H and J are therefore of iden¬ 
tical form in a good conductor. We now turn our attention to solving 
for E: the other solutions follow the same pattern. 
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17.3 Field penetration into a good conductor 

• The following analysis will show that the E-field (and, similarly H 
and J) only penetrates with significant amplitude to a limited depth 
into a good conductor. At high frequencies, it will be found that this 
depth is small in comparison to the wavelength. This implies that the 
performance of high frequency fields near a conductor is not affected 
much by the shape, edges or curvature of the conductor, because typical 
conductor dimensions are a lot higher than the penetration depth of 
the field. We therefore study an E-field of magnitude E 0 at the surface 
of a conductor, and can make the assumption that the conductor is 



in the ^-direction and which propagates in the x-direction. There is 
therefore no variation of electric field with y or z. As this means that 
= fz = 0, our differential equation for E, namely Y 2 E = jufiaE, 
becomes 

cP E ~ ~ /v 

= jui[xcrE z = t 2 E z (17.15) 

for which 

r = (17.16) 

where the factor is recognised as the positive square root of j, / is 
the frequencv in hertz, and 6 = , 1 in metres. 
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• The differential equation for E z can now be written 

d 2 E z , - 

~d^~ TEz = ° (17J7) 

and its characteristic equation, s 2 — r 2 = 0, implies that 

(5 + r)(s — r) = 0 and so s = ±r. The solution is therefore 

E z = Ae~ TX + Be TX (17.18) 

where A and B are constants to be determined from the boundary 
conditions. We quickly see that, since E z certainly cannot increase 
to infinity as x increases, B must equal zero and since, when x = 0, 
E z = E 0 , A must equal E 0 . Hence our solution for the differential 
equation is the decaying exponential 

E z = E 0 e~ TX = E 0 t~^) x = E 0 e-*‘e- j f (17.19) 

• Likewise, you could arrive at H y = H 0 e~?e~*% or J z 
very similar developments. 

• The interpretation of the E-field solution is shown 
in this figure. You will see that the magnitude of the 
z-component of E is E 0 /e at depth x = 6, which is 
known as the skin depth, and that the envelope of 
the magnitude of the E-field decays exponentially 
with depth. Furthermore, the decaying sinusoidal 
wave associated with the E-field travels into the 
conductor in the x-direction. 

• Note also that, in all cases the E (or D) and H (or B) 
tion which is uniquely determined by the right-hand set [E, H, x], where 
x is in the direction of propagation of the wave, and E = [E x , E y , E Z ] T . 


-J 0 t < e }X 6 by 



fields have direc- 
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17.4 Resistance and internal impedance 

• Let us now consider a section of the semi-infinite plane of unit width and 
length (but effectively infinite depth). We can calculate the total cur¬ 
rent flowing in this conductor by integrating the exponentially-decaying 
current density down to infinite depth. Thus 

Itot = f°° J 0 e~^ x dx = (17.20) 

Jo 1 +J 

• The magnitude of the surface electric field is then simply found from 
the surface current density to be: 

E 0 = — (17.21) 

<7 

The internal impedance of the conductor per unit length and width is 
the ratio of the field and of the current at the surface, which is 

Z s = R s + jto Li = — = — + (17.22) 

Itot &a 

• - Note that this expression gives us the surface resistivity , R s , of the 

conductor ( multiply it by length and divide by width to get the 
resistance of the conductor), and also the internal reactance , 

The latter does not include any contribution to total inductance 
from fields that are external to the conductor. 

— Note also that for a plane conductor of unit length and width, 
the resistance and internal reactance are of equal magnitude or, 
equivalently, that the phase angle of the internal impedance, Z s , 
is 45°. 

— Finally, note that the resistance of the conductor is inversely pro¬ 
portional to skin depth, 8, which is itself inversely proportional to 
the square root of the frequency, /. Thus, at high frequencies we 
understand that fields penetrate less deeply into a conductor, and 
that the resistance of the conductor increases. This is in line with 
our expectation for a conductor of decreased cross-sectional area. 






Chapter 18 

The Basis of Kirchhoff’s Laws 

18.1 Lumped elements 


• In everyday electrical engineering, the great majority of electromagnetic 
interactions are analysed and evaluated using standard circuit theory. 
The purpose of this chapter is to show that all circuit theory has its 
foundations in field theory, and that therefore the use of circuit analysis 
techniques is justified upon solid theoretical grounds. 

• In lumped-element circuit analysis, we make the assumption that all 
of the circuit elements used are small compared with the wavelength 
of the electromagnetic fields within them. This assumption allows us 
to treat the conditions around each element as quasistatic , simplifying 
our treatment considerably. At high frequencies (e.g. in microwave 
circuits) or where very large circuit elements are being considered (e.g. 
an inter-city power transmission line), the quasistatic assumption is 
invalid. In such cases, it is necessary to resort to transmission-line 
theory in preference to standard circuit analysis. We concern ourselves 
here only with the validity of the latter at sufficiently low frequencies. 

• The three ideal lumped elements most frequently dealt with by circuit 
analysis (see Module A of EEE221W) are the resistor , the capacitor 
and the inductor. We now consider each of these in turn from the 
standpoint of field theory. 


Ill 
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18.2 Resistors 

T 

'*- » . ? 

j7 

• The resistor is a dissipative element: that is, the interactions which 
occur within it at a subatomic level cause it always to absorb net elec¬ 
trical power from its surroundings. It converts much of the energy that 
it absorbs to the generally unusable form of heat. 

• Resistors conduct current, and are characterised by a conductivity , cr 
in Um -1 . Hence, as is well-known, they obey Ohm’s Law, so that if 
any field E is applied across a resistor, then a current of density J will 
flow, where 

J = crE ( 18 . 1 ) 

• The voltage across the resistor can be expressed using the simple rela¬ 
tion for electrostatic potential difference: 

f 2 f 2 J f 2 I dl 

V 2 i = $2 - $1 = - / E • dl = - / - • dl = - / — ( 18 . 2 ) 

J 1 Ji a J\ a A 

where it is clear that J || dl, that J = I/A, and that A represents the 
cross-sectional area of the resistor. 


• If we now assume that the current, 7, the conductivity, a and the 
cross-sectional area, A, are all constant between points 1 and 2, we can 
write 

V 21 = —[ 2 dl= ~— a = -IR (18.3) 

<7 A J\ cr A 

where L is the length of the resistor, and R = L/(a A) is known as the 
resistance. We have thus recovered the familiar form of Ohm’s Law. 


• In the case where 7, cr and/or A are not constant over the length of the 
resistor, we simply define 



2 dl 
cr A 


(18,4) 


to arrive at 


V21 (t) = —Ri(t) 


(18.5) 
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18.3 Inductors 



J7 

>■ 

i 

9 

✓ 

/ 


• Inductors, as we know, store energy in the magnetic field that is asso¬ 
ciated with the current that passes through them. Ideal inductors have 
zero resistance in their coils, and so all of the stored energy may be 
recovered at a later time. 


• Taken around the path shown in the diagram, clockwise from point 2, 
the integral of the electric field around the path is evidently 


j E ■ i\ = F E ■ i\ + I E • il 


(18.6) 


where the first term on the right-hand side integrates the field along 
the perfectly-conducting coil, and is therefore equal to zero, while the 
second term integrates the field along the return path between the 
terminals. We thus have, by Faraday’s Law, 

2 o 

j>E ■ d\ = 0 + E ■ dl = —— ■ dS (18.7) 


• We defined inductance much earlier in the course as 

L = \ [ B • dS 
I Js 

so we may now write, for the quasistatic case, 


(18.8) 


**<'> = - r E - 1,1 = I /, B • <« - (18 - 9 > 


if the current through the inductor is taken to be the time-varying z(t). 
Again, this is a very familiar result from circuit theory. 
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• Note that if the resistance of the coil were not zero, but instead had 
value Rl , we would then have 

fE-d\ = +R L i(t)-v 23 (t) = -^J s B-dS (18.10) 

so that, as we might expect, a time-varying current would now give us 

di(t) 


v 23{t) — R L i(t) + L- 




dt 


(18.11) 


V xr 


Y 


• Note also that a voltage may appear across an inductive coil in a circuit 
as a result of magnetic flux changes that are due to some other source 
of magnetic flux. Thus, if coil Y lies close to coil X and has a time- 
varying current iy(t) flowing in it, then the magnetic flux that links 
coil X due to this current is written ipXY- By direct analogy with the 
concept of self-inductance, L , the mutual inductance , M, of the coils is 
defined as the ratio of the flux linkage and the current in coil Y that 
produces it. Thus 

M XY = ^ (18.12) 

ly 

The voltage induced across coil X by this flux linkage is then written 

(18.13) 




dt 


dt 


The phenomenon of mutual inductance is symmetrical in the sense that 
the same value for M is used in calculating the voltage in coil Y due 
to changes in the current in coil X. Note that dots drawn at the ends 
of the coils in the circuit diagram show the directions of the windings, 
and so give the polarity sign for the voltage induced. 
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18.4 


Capacitors 

3 

#- 


Q 


X o 

■» » 


• Capacitors store energy in the electric fields between their plates, and 
we established much earlier on that, with a charge Q on either plate, the 
capacitance is found using the potential difference between the plates, 
as 

c _ Q _ Q 

$3 ~ $0 u 3o(^) 

• Any current flowing in the capacitor will, by continuity, merely raise 
or lower the value of Q. Thus we can write 

!(,) = f = = cdJ 7T < 18 ' 15 ) 

or. equivalently, 

v 3 o{t) = 77 / i(t)dt (18.16) 

C Jo 

if the capacitance is constant over time, and if the capacitor is initially 
uncharged. We have therefore recovered the familiar expressions for 
the behaviour of an ideal capacitor in a low-frequency circuit. 

18.5 Kirchhoff’s Voltage Law 


• We can now put together the three elements considered in this chapter, 
and drive them with a time-varying voltage source, as shown in the 
diagram below. 
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Taking the line integral of electric field around the whole path, clock¬ 
wise from point 0, now gives us 

<f E • d\ = /' E • d\ + [ 2 E • d\ + / 3 E • d\ + /° E ■ d\ = / B • dS 

J Jo Ji J 2 J 3 ot Js 

(18.17) 

• The surface 5 referred to in the equation above is, of course, any surface 
bounded by the loop of the circuit. If we assume that the magnetic 
flux in the inductor is largely confined to the inside of the coil (a fair 
assumption, as we know from our study of the solenoid), we can see 
that the component of stray magnetic flux perpendicular to the surface 
S is negligible. Insofar as it needs to be considered at all, it could in 
any case be added to the effect of the inductor, L. Hence, the term on 
the right-hand side of the above equation can be regarded as equalling 
zero. The equation thus becomes 

— j) E • dl = Vio + t’ 2 i(*) + v 32 (t) + ^ 03(0 (18.18) 

= v 0 {t) - Ri(t) - - 77 t i(t) dt (18.19) 

di 0 J 0 

= 0 (18.20) 

• This is what we know as Kirchhoff’s Voltage Law, which we use habit¬ 
ually in electrical circuit analysis. You can now do so with confidence, 
knowing that it is a result which is formally derivable from electromag¬ 
netic field theory. The general statement of the law for N elements 
around a closed path in a circuit is 

N 

n=l 


(18.21) 
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18.6 Kirchhoff’s Current Law 



• Consider a node in an ideal circuit at which a number of conductors 
(assumed perfect) all meet. To be completely general, let us suppose 
that there is a charge Q at the node: evidently, this charge must be 
maintained by a capacitive effect, or else it would instantly be dissi¬ 
pated in the form of currents along the perfect conductors. 

• We know that the principle of continuity states that 

V J = -| OT i 3dS = -jtl<’ dV < 1822 > 

and we can give meaning to these expressions in this situation by imag¬ 
ining setting a surface 5 around the node, enclosing volume V, which 
would then have charge density p = Q/V as shown above. In the 
integral form of the continuity expression, 

i 3dS+ mi l,dV = 0 (18 ' 23) 

the first term clearly represents the sum of the currents, /), / 2 , h, I 4 ... 
which flow in the perfect conductors. The second term is dQ/dt, the 
current leaving the node via the capacitor. This current is not included 
in the first term, since J denotes conduction or convection current only, 
and the capacitor current is a displacement current. 

• Nevertheless, the equation shows that the sum of all the currents of all 
types out of the node is zero. We know this as Kirchhoff’s Current Law , 
and we use it implicitly on a daily basis when we work with circuits. 
In general form, it may be written 

N 

X> = ° 

n —1 


(18.24) 
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• With Kirchhoff’s Laws established from field theory, you may recall (see 
EEE221W) that the whole structure of circuit theory now follows from 
them. We therefore have the theoretical basis upon which to analyse 
circuits and their equivalents, thanks to the theory of electromagnetic 
fields. To be specific, Kirchhoff’s two laws are necessary to the proofs 
of the following well-known circuit analysis techniques: 

— Nodal analysis 

— Mesh analysis 

- The theorems of Thevenin and Norton 

- AC phasor analysis 


Chapter 19 

Reflection of Travelling Plane 
Waves 

19.1 Intrinsic impedance of a medium 


• In introducing the concept of plane electromagnetic waves in free space, 
we saw how Maxwell’s equations reduced (with J = p = 0) to 

VxE = -/i- and V x H = e— (19.1) 

dt at 

We can extract component expressions from these equations, and, for 
a plane EM wave propagating in the z-direction, these are: 


dEy _ 

dH x 

dHy 

dE x 

dz 

= » dt 

dz 

~ _e dt 

dE x 

dHy 

dH x 

dEy 

dz 

= “at 

dz 

~ e dt 

0 

dH, 

0 = 

dE z 

= dt 

6 dt 


(19.2) 

(19.3) 

(19.4) 


Further study of these equations quickly led us to the conclusion that 
the E and H fields of a travelling wave are transverse to the z-direction 
in which the wave propagates, and also led us to the wave equation, 
which we are now able to solve for any of E x , E y , H x or H y . 
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• We then found, for example, that the general solution to the wave equa¬ 
tion for E x gave us a wave travelling in the direction of +z (which we 
will now denote E x+ ), and also a wave of the same velocity propagating 
along —z (which we can similarly call E x -). Thus the solution was 

E x = E x+ + E x _ = /i (^t — — ^ + f 2 ^ t + —^ (19.5) 

where v = 1 /y'/Ie was the velocity of the wave, which equals the vacuum 
speed of light, c, in free space. 

• In the commonly-encountered situation where the time-variation of the 
fields is sinusoidal (with frequency u>), we used phasors and solved the 
Helmholtz equation to arrive at the similar solution for E x : 

E x = E x+ + E x - = Cl e~ jkz + c 2 e jkz (19.6) 


where 


iO 


k = iOy/JIe = — 


• Now, one of our component equations above states that 


dE x 

dz 


= -I*- 


dH v 


dt 


which, in phasor form, becomes 

dEx 

dz 


so that 


= -njuHy 


H y = 


1 dE x 
ju/i dz 


L0/J, 


Cl e~ jkz - c 2 C kz 


E x + E x — 


(19.7) 

(19.8) 

(19.9) 

(19.10) 

(19.11) 

(19.12) 
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• Hence we can write H y as 

Hy = H y+ + Hy- = ^±-^ (19.13) 


where rj = y/(fj,/e) is therefore the ratio of E x (the x-component of E- 
field) and H y (the y-component of H-field) everywhere along the wave. 
You can see that p is a constant for any given medium, and, because 
it has dimensions of (V/m)/(A/m)= fl, it is known as the intrinsic 
impedance of the medium in which the wave propagates. 

• In particular, the impedance of free space is 

= 376.73 « 120tt fl (19.14) 

This is a quantity with which you will frequently meet in any analysis 
of waves travelling in space or in similar media. 

• As we remarked when studying the wave equation, it is a simple pro¬ 
cedure to extract from Maxwell’s equations the wave equation for any 
of the other quantities E y , H x or H y . You should practise doing this. 
You will then find that, if you have the solution for E y such that 

E y = E y+ + E y . = c z e- ]kz + c 4 e :kz (19.15) 



then a resubstitution into a component equation, similar to the one 
described above for E x , will give you 


H x = H x+ + H x . = -■ 


J y+ 


+ 


E v - 


V 


V 


(19.16) 
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• Combining the results of both procedures, we can now see that 


Ex+ 

Hy+ 


E'y± 

H x+ 


= V 


and 




(19.17) 


One important interpretation of these equations is that they show that 
the total E and H fields of a travelling plane wave are everywhere per¬ 
pendicular to one another. This is illustrated in the following diagram, 
which depicts an electromagnetic wave moving out of the page (in the 
positive z-direction). Once the information contained in equation 19.17 
has been included in the diagram, you can see by a similar-triangles 
approach that E and H must be perpendicular to each other, and also 
that the wave travels in the direction of E x H, as expected. 
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19.2 Incidence on a perfect conductor 


Sp 




ro-Hec 


We now imagine a plane wave, that has been 
propagating in a medium such as free space, - 

impinging upon a perfectly-conducting surface. 

In the simplest case, we can visualise the wave 
striking the conductor so that its direction of 
propagation is normal to the conductor’s 
surface (i.e. it is in the z-direction as drawn). 

Now, at the surface of the perfect conductor 
there can be no x- or y- components of E-field, 
because such a field would instantly collapse as 
charge carriers moved in response to it. Hence, the negatively-travelling 
part of the plane wave must, at the surface, possess x- and y- compo¬ 
nents of E-field which exactly balance (i.e. they are equal and opposite 
to) the components of the impinging, positively-travelling part of the 
wave. Under these circumstances, we would refer to the negatively- 
travelling part of the wave as a reflection ; and, since its components 
are equal in magnitude but opposite in direction to those of the incident 
wave, it is clear that the reflected wave is very similar to the incident 
wave, but that it moves away from the conducting surface, as expected. 


-»Z 


Next imagine a plane wave impinging upon the 
conductor at some angle 6 to the normal line 
as shown here. Again, the x- and y- components 
of E-field at the conductor’s surface must 
each total zero, and so we expect there to be a 
reflected travelling plane wave, which is drawn 
here at an angle of O' to the normal. 



«♦ 


In this example, we have drawn the H-field 
coming out of the paper, and the E-field lies in the 
plane of incidence , which is the plane containing both 
the direction of propagation and the normal line (it is 
the plane of the paper as drawn here). This geometry is 
referred to as transverse magnetic , or parallel, polarisation. 



J 
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Note how the diagram preserves the perpendicularity of E and H, and 
note also the fact that E x H gives the direction of propagation for 
both the incident and the reflected waves. 

• Now, the directions of the incident and reflected waves have been de¬ 
noted ( and (—)£' respectively, so we can write the whole E-field of the 
wave as . 

E = E+e~ jk< + E-e jk( ' (19.18) 


The diagram also shows how we can write ( and (' in terms of x and z 
as follows: 

( = x sin 6 + z cos 9 and (' = — x sin O' + z cos 6 ' (19.19) 

Hence, the whole E-field of the wave may be given as 

E (x,z) = E + e~ 3k(xsin6+zcose) + E_e M ~ xsine ' +zcos9,) (19.20) 


Just considering the x-component of E-field, and still with reference to 
the diagram, you can easily see that taking x-components gives 


E x (x,z) = -E + C os0e- J * (xsin0+2COS(?) + £_ cos 0 ’ e M-*s™o'+zcose') 

(19.21) 


• At the boundary with the perfect conductor z = 0, and we know that 
E x must equal zero as well. So, setting E x (x, 0) = 0 gives us 


E + 


cos 0 e- 3kxsine = £_ cos 6 'e- jkxsia0. 


(19.22) 
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The only wav that the complex numbers on both sides of this equality 
can be equal for all values of x (i.e. anywhere along the surface of the 
conductor) is if they are identical in both magnitude and phase. We 
conclude, therefore, that 

k sin 6 = k sin 9 so 6 = 8' and E + = E- (19.23) 


and so we see that the reflection at the perfect conducting surface is 
lossless, and that the angle of reflection to the normal is equal to the 
angle of incidence. 


• You should not find it difficult to come to the same conclusion in study¬ 
ing this second diagram, in which the incident and reflected waves both 
have their E-fields directed straight into the paper (perpendicular to 
the plane of incidence). This is referred to as transverse electric , or per¬ 
pendicular, polarisation. Since any other orientation of the E-field may 
be regarded simply as a superposition of parallel and perpendicular 
polarisations, we know that the important conclusions that reflection 
is symmetric and lossless will be true, irrespective of what linear po¬ 
larisation the incident wave has. 
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19.3 Incidence on a dielectric boundary 

• Unlike at the boundary of a perfect conductor, across which a plane 
wave cannot transmit energy, a travelling wave is not so constrained at 
the boundary between two dielectric media. At such a boundary, we 
know from our early consideration of boundary conditions that the tan¬ 
gential components of the E and H fields are unchanged by the change 
of medium. Furthermore, there is some transfer of energy from the first 
medium to the second as the wave crosses the boundary; and yet not all 
of the incident energy is transmitted because the normal components 
of E and H field are not the same in both dielectrics. Hence, we expect 
to see both transmission (known as refraction) and also reflection in 
such a situation. 

• Following our approach at the perfect conductor 
boundary, and with reference to the diagram, we can 
write an equation for the x-components of E-field at 
a point (x,0) anywhere on the dielectric boundary: 

? 

E x (x, 0) = -E+ cos 0 ie -^* sin *i +£_ cos e\t~ 3k ' xsin<) ' (19.24) 

= -E 2 cos6 2 e- jk2Xsin92 (19.25) 

Since the boundary condition dictates that 

E n = E t2 (19.26) 

arid, since the previous expression must be true irrespective of the value 
of x (i.e. everywhere on the boundary), we can conclude, similarly to 
the case of the perfect conductor, that the following terms, known as 
phase factors, must be equal: 




(19.27) 


k\ sin 9\ = k\ sin = k 2 sin 0 2 
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• We can obtain two very useful facts from these equations. The first 
equality implies, just as with the conducting boundary, that 


0i = 0J (19.28) 


so that reflection is symmetric about the normal line, much as we might 
expect (consider air and water as the two dielectrics, for instance). The 
second equality, together with several of our earlier definitions, yields 


ni sin 02 k\ 
n 2 sin 0i k 2 



(19.29) 


where 



(19.30) 


is called the refractive index of medium x, and the approximation signs 
hold because p x rs ~ Ho for all common media. You should recog¬ 
nise that we have therefore proven Snell’s Law of Refraction from the 
standpoint of field theory. 


• We can also obtain useful expressions for how much of the incident wave 
is reflected and how much is refracted at the dielectric boundary. To do 
this, we define a quantity called the wave impedance , which is simply 
the ratio of the electric and magnetic fields in particular directions. 
Thus, in the first medium, the wave impedance is 

Z\ = 7W" “ - TT~ = 0 1 cos 0 J (19.31) 

tt y + tty- 

Likewise, on the other side of the boundary, we imagine that the wave 
“sees' 1 a load of impedance 




Ex 2 

Hy 2 


(19.32) 
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• Now, the continuity of the tangential components of E and H at the 
boundary gives us 


E x + + E x - — E x 2 and H y + 4- H y ~ — H y 2 (19.33) 


It should be clear that the second of these two equations gives us the 
equation below. (We will then substitute the first equation into this new 
equation in two different ways to quantify the reflection and refraction 
at the boundary). 

E x + E x - E x 2 



1. Using the first equation to substitute for E x2 gives us 


so 


E x + E x - 


E x + E x _ 

zT zT 


(19.35) 


E x + 


■ 1 1 ■ 

= E x . 

• 1 1 ' 

■ Zj - Z L . 




(19.36) 


and from this we can define the reflection coefficient , p, which 
turns out to be 


P = 


E x - 

E x + 


j_i_ 

Zl_ Zl 


-L + J- 

^ Z L 


Zl-Z x 
Zl + Zi 


(19.37) 


Hence, if we know the wave impedances of the two media, we 
can predict what fraction of the incident wave is reflected at the 
boundary. 
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2 . In a similar way, using the first equation to substitute for E x - will 
give 


E x + ( E x 2 E x +) E x 2 

Z\ Z\ Zi 


(19.38) 


which quickly gives 


2 E x + 
~Z^ 


= E. 


x2 


-Z[ + Z~i. 


(19.39) 


so we can then define the transmission coefficient , r, as 


E X 2 _ Zi _ 

-£*+ zL + ZI Z\ + Zl 


(19.40) 


This can be used to predict how much of the incident wave is 
refracted at the boundary. 


• You will certainly have noticed the similarity between these equations 
and those which you have met in your study of transmission lines. This 
should confirm for you that the guided transmission of electromagnetic 
waves that one studies in transmission line theory is described by field 
theory just as effectively as the case in which the waves propagate, 
unguided, through the surrounding medium. 



- 
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Test 1 Specimen Questions 


1995 Test 1 


Attempt any three questions, and state all your final answers 
in engineering units. Credit is given for showing all of your 
working carefully. The time allowed is 50 minutes. 

60 marks will gain 100% 

1995 Test,! Questionl (20 MARKS) 


(a) Show that the electric field intensity due to a ring of positive charge, of 
charge density pi and radius a, perpendicular to the z-axis, is given at any 
point on the z-axis with co-ordinate z by 


E = 


Pi az 

2e(a 2 + z 2 )t ez 


[ 8 ] 


(b) Hence state the magnitude and the direction of the electric flux density 
vector, D, at a distance of z = 10cm up the axis from the centre of a charge 
ring of radius 10cm, assuming a charge density of pi = 2Cm -1 . Be careful to 
state your answer with units. [6] 


(c) Make a neat sketch of the electric field in the vicinity of the charge ring, 
showing the behaviour of the field lines both on and off the ring’s axis. What 
is the appearance of the field at a distance from the ring that is much greater 
than its radius? [6] 
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1995 Testl Question2 (20 MARKS) 


(a) Calculate the electric field intensity between the two cylinders of a coax¬ 
ial line separated by a medium of permittivity t. Use Gauss’s Law, with a 
special gaussian cylindrical surface of length l and radius r to do this. You 
should take the charge per unit length on the inner cylinder as /);, and should 
obtain the answer 


E r 


_Pl_ 

2tt er 


[5] 


(b) Now calculate V • D in the region between the cylinders. Is the result 
what you might have expected? Explain. [5] 


(c) The inner cylinder is at a potential of IV and has a radius of 1.831cm. 
The outer cylinder is at a potential of 5V and has a radius of 4.978cm. Use 
this information to solve Laplace’s equation between the two cylinders, to 
arrive at the following formula for the distribution of the potential at any 
radius r: 


$ = 4 In |r| + 17 


[5] 


(d) From the answer to part (c), find the radii of the circles of equipotential 
at 2V, 3V and 4V. These circles have been drawn to scale for you on the 
accompanying full-scale diagram. Add equally-spaced field lines to this dia¬ 
gram, and hence deduce the capacitance per metre of the coaxial line. 


[5] 
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1995 Testl Question3 (20 MARKS) 


(a) For two point charges q and — q at [d/2, 0. 0) and [—d/2, 0,0) respectively, 

find the potential for any point [x,y, z) and from this derive the electric field. 
Check the results by adding vectorially the electric field from the individual 
charges. [4,4,2] 

(b) Find the electric field and charge density as functions of x for a space- 
charge-limited, parallel-plane diode with potential variation given by 

* = 1 5 ' 5 ) 

1995 Testl Question4 (15 MARKS) 


(a) Show that a parallel-plate capacitor with plates of area A and separation 

d, charged to and at a charge density of ps Cm -2 , has a flux density 
between the plates of magnitude D = ps . State clearly any assumptions that 
are made in your explanation. [4] 

(b) Use the definitions of potential and capacitance to show: 

(i) that the voltage between the plates is V c — ps d/t\ [5] 

(ii) that the capacitance of the parallel plates is C = tA/d [5] 

(c) Then, given the result that the energy stored in an infinitesimal volume 
dV of an electric field E of flux density D is 

dU c = • EdV, 

show that the energy stored in the whole field of the capacitor is 

£4 = \cv 2 

Again, state clearly any assumptions made in your analysis. 


[ 6 ] 
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1996 Test 1 


Attempt any two or three questions, and state all your final 
answers in engineering units. Credit is given for showing all 
of your working carefully. The time allowed is 50 minutes. 

50 marks will gain 100% 


1996 Test! Questionl (20 MARKS) 


(a) Show by integration of the contributions from each differential charge 
that the electric field intensity due to an infinite straight line of positive 
charges, of uniform charge density pi distributed along the z-axis, is given in 
cylindrical co-ordinates at a distance r from the line by: 


E 


Pi 


2~Ktr 


where t is the permittivity of the surrounding medium and r is the unit radial 
vector. [7] 


(b) Verify the result given above by applying Gauss’s Law to a carefully 
chosen gaussian surface near the line of charges. Present your argument in 
terms of a clearly-labelled, neat diagram, and pay particular attention to 
correct notation. Explain your steps in words where necessary. [7] 


(c) Although the result in part (a) above refers to a line of static charge, it is 
also adequate for finding the electric field due to a steady current in a wire. 

(i) Explain briefly why this is so. [3] 

(ii) Make a sketch of the electric field near a wire carrving a steady current. 

[3] 
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1996 Testl Question2 (20 MARKS) 


A parallel-plate capacitor consists of a pair of large conducting plates which 
are separated by a dielectric medium of relative permittivity e r = 3.0. The 
plates are perpendicular to the 2 -axis, as shown in the upper drawing on 
page 6, and each plate is a square of side 5cm (i.e. of area 0.0025m 2 ). The 
plate separation is d = 5mm. 

At a given instant in time, currents in the wires leading to the capacitor have 
produced an accumulation of positive charge on the upper plate and negative 
charge on the lower plate, so that the plates have charge density p s Cm -2 
and — p s Cm -2 respectively. 

(a) Given that the electric field intensity due to each plate has magnitude 

E = briefly state (with a reason) the total electric field intensity vector in 
the region between the two plates. [2] 

(b) From your answer to part (a) above, show that the potential function 

between the capacitor plates, as a function of the distance from the upper 
plate, is $ = 6 — using the fact that at the upper plate (where 2 = 0) 
the potential is 6V. [ 4 ] 

(c) Given also that the lower plate is at ground potential (0V), use the result 
given in part (b) above to state the total charge Q on the upper plate. [ 4 ] 


(d) A cross-sectional view of the parallel-plate capacitor is drawn to scale on 
page 6, (which you must hand in, if you answer this question): 

(i) Add to this drawing the equipotential lines for $ = 2V and $ = 4V, 
and also add appropriate field lines to complete a graphical field map. [3] 

(ii) Calculate the total capacitance of the plates from your field map, 

ignoring the effects of fringing (i.e. consider only the field lines and the 
equipotentials that are directly between the plates). [3] 


(iii) Compare your answer to part (ii) above with the capacitances as cal¬ 
culated from the well-known expressions: 


C = 


A 



and 


[ 2 , 2 ] 
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Question 2(d) NAME: 

(i) 



(ii) 


(iii) 
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1996 Test! Question3 (20 MARKS) 

A steady current in a long wire gives rise to an electric field E = 
where pi is the net charge density of the current carriers which are uniformly 
distributed along the wire, e is the permittivity of the surrounding medium, 
and r is the radial distance, in cylindrical co-ordinates, from the central pr¬ 
axis of the wire. 

(a) Show that the divergence of the electric flux density vector is equal to 

the net charge density both inside and outside the wire. [3,3] 

(b) The wire can be viewed as an ohmic resistor of conductivity a. State: 

(i) Ohm’s Law, in terms of field quantities; 

(ii) The Principle of Continuity in differential form; 

and use these equations to show that the potential function $ in the vicinity 
of the wire will obey Laplace’s Equation: 

V 2 $ = 0 [2,2,2] 

(c) In the solution of Laplace’s Equation in this case, explain very briefly 
why it is acceptable to neglect: 

(i) the tangential {<j>) component of the potential; 

(ii) the longitudinal ( 2 ) component of the potential. [ 1 , 1 ] 

(d) Solve Laplace’s Equation for this case, with the boundary condition that 
at a distance r 0 from the wire the potential is zero, and hence show that the 
potential function around the wire as a function of the radial distance r is 

$ = A In |r| — A In |r c | 

where A is a constant. [4] 

(e) Calculate the potential function around the wire directly from the given 
expression for the electric field intensity, and hence show that the constant 
A in part (d) above is given by 
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1996 Testl Question4 (20 MARKS) 


(a) Write down the relationship between electric field intensity and electro¬ 
static potential both in its integral form and in its differential form. Start 
from this relationship to explain in careful detail why the equipotentials are 
everywhere perpendicular to the field lines in an electrostatic field. [ 1 , 1 , 3 ] 


(b) Write down Gauss’s Law both in its integral and in its differential form, 
and show that these two expressions of Gauss’s Law are equivalent. [ 1 , 1 , 3 ] 


(c) Give a simple verbal statement of what is meant by the divergence of a 
vector quantity, and then derive from first principles the expression for V • D 
in the cartesian co-ordinate system, namely: 


V D 


dD x dD y dD z 
dx dy dz 


[ 8 ] 


(d) A research engineer calculates the electrostatic potential function at a 
moment in time between the stator and the rotor of a motor. The function 
is very complicated, and the researcher is obliged to find it by using an ex¬ 
pensive computer package which applies elegant numerical methods to the 
boundary conditions. These boundary conditions have to be input for the 
particular case of the motor in question. Once the computer completes the 
calculation, what does the engineer know for certain about the solution for 
the potential function that the computer gives? 
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1995 Test 2 

Attempt any two questions in full plus part (a) only of the 
other two questions. State all final answers in engineering 
units. Credit is given for showing all of your working carefully. 
The time allowed is 50 minutes. 48 marks will gain 100% 


1995 Test2 Quest,ionl (20 MARKS) 


(a) State the Biot-Savart Law, and provide a diagram which illustrates all of 
the seven variables that are connected by this Law. [4] 


(b) Use the Biot-Savart Law to calculate the magnetic field vector (magnitude 
and direction) at a point distance r from a straight wire of infinite length 
which carries a dc current 1 0 . You may take as given the integral 



2 

r 2 ^/(r/a) 2 + 1 


[5] 


(c) State Ampere’s Circuital Law, and explain carefully (5-10 lines plus a 
diagram) how this Law yields the same result as your answer to part (b) for 
the magnetic field of an infinite wire carrying dc current I 0 . [5] 

(d) From the defining equation for magnetic flux density, di — ld\ x B, use 
the result of parts (b) and (c) to calculate the magnitude and direction of the 
force per metre length between two parallel wires lm apart, each carrying a 
current of 1A in the same direction. (This is often regarded as the primary 
definition of the ampere in the SI system). [6] 
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1995 Test2 Question2 (20 MARKS) 


(a) Make a sketch of the magnetic fields lines in the vicinity of a parallel- 
microstrip transmission line. [4] 


(b) A microstrip transmission line consists of two long parallel plates of width 
w, copper thickness h and separation d. The plates are separated by a 
dielectric of magnetic permeability fi. Listing three important and distinct 
assumptions that are made in your analysis, use Ampere’s Circuital Law to 
show that the inductance per unit length of such a transmission line is 

L =^ [Hm -1 ] [7] 

W 


(c) With reference to your sketch in part (a), explain what you would expect 
to happen to an imaginary “curl meter paddle” placed exactly between the 
two conducting plates of the transmission line discussed in part (b). Can you 
justify this expectation in terms of a mathematical expression for the curl of 
magnetostatic field? [5] 


(d) Show that, if the magnetic field between the parallel plates is written H 0 
and current 1 flows in them, then, given the assumptions listed in part (b), 

JA h ^ v =\ li2 


[4] 
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1995 Test2 Question3 (20 MARKS) 


(a) What is meant by describing a magnetic material as 

(i) linear (ii) isotropic 

It is acceptable to give each answer either as a simple equation or as a very 
short explanation. In either case, give your answer in terms of the quantities 
B, H and p. [4] 

(b) Starting with the expression for stored magnetic energy, 

dU H = ! H dB dV 
Jv 

explain how and why it is possible to estimate the energy dissipated in a 
unit volume of a non-linear , isotropic , hysteretic material, given a plot of the 
BH response of the material over a complete cycle of periodic magnetic field 
variation. [8] 


(c) The diagram below is an approximation of the hysteretic loop for fairly 
pure iron. Use the result that you gave in part (b) to show that, in a sample 
of this material of volume 27r x 10~ 3 m 3 , 5J of energy is dissipated per cycle 
of the periodic magnetic field variation. You are given that 


• 1 gauss= O.lmT 

• 1 oersted= 250/7rAm -1 

• 1 tesla— lNm _1 A _1 



[ 8 ] 


(4,ttooo) 


H (otrst^Js) 
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1995 Test.2 Question4 (20 MARKS) 


(a) Write down Maxwell’s four equations for time-varying electromagnetic 
fields, in their differential form. [4] 


(b) EITHER: Explain in detail how Faraday’s Law of Electromagnetic In¬ 
duction is used to obtain the expression for the curl of electric field in the 
time-varying case, as a generalisation of the static case in which V x E = 0; 

OR: Explain in detail how the concepts of continuity and of displacement 
current are used to obtain the expression for the curl of magnetic field in the 
time-varving case, as a generalisation of the static case in which V x H = J. 

[ 8 ] 


(c) Show that your four answers to part (a) are equivalent to the integral 
forms: 

1. f s D-dS - f v p dV 

2. f s B-dS = 0 

3. / E • dl = — Jj / s B • dS 

4. f H • dl = / 5 J • dS + / 5 D • dS 

and state very briefly what each equation expresses, in either its differential 
or its integral form. [8] 
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1996 Test 2 


Attempt any three questions. Credit is given for showing all 
of your working carefully. The time allowed is 50 minutes. 

45 marks will gain 100% 


1996 Test2 Questionl (15 MARKS) 

On part of a certain printed circuit board, a dc current of 500mA flows in 
a semicircular loop of radius a = 4cm, as illustrated below. Show that the 



point with distance 2 from the radial centre of the semicircle along the 2 -axis 
of the semicircle (this axis being perpendicular to the printed circuit board), 
has c-component 

la 2 


H x = 


and y-component 


Hy = 


4(a 2 -f z 2 )'? 


laz 

2ir(a 2 + 2 2 ) 2 


[4] 


and x-component 


[4] 


H x = 0 


[ 1 ] 

Hence state the magnitude and the angle to the 2 -axis of the magnetic field 
intensity vector at a point with distance 2 = 3cm along the 2 -axis from the 
semicircle’s radial centre. 
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1996 Test2 Question2 (15 MARKS) 

It is common practice in some industrial applications to suppress unwanted 
signal spikes in a conductor by threading the conductor through a cylindrical 
bead as shown below. 



(a) Explain carefully how you would calculate the inductance due to the 
bead in this configuration, and derive a formula for this inductance 
in terms of the bead’s dimensions, a, b and /, and in terms of the 
permeability, /u, of the material of which the bead is made. 

[5] 

(b) Use the formula which you derived in part (a) to calculate the induc¬ 
tance due to an actual bead which has an inner radius of 1mm, an outer 
radius of 2.15mm, a length of 7.2mm, and which is fabricated from a 
material of relative permeability = 900. 

[4] 

(c) The actual inductance of this bead has been measured accurately, and 
was found to be 1//H. Can you explain any small discrepancy between 
this measured value and your calculated answer above? 

[3] 

(d) Describe how and why the bead acts as required in the suppression of 
unwanted signal spikes. 


[ 3 ] 
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1996 Test2 Question3 (15 MARKS) 


(a) Show clearly that electrostatic fields have zero curl. 

[3] 

(b) Starting from the basic definition of curl, use Ampere’s circuital law to 
derive an expression for the curl of magnetic field intensity. 

[5] 

(c) Magnetic flux density may also be found from the curl of magnetic 
vector potential. Use this fact together with Stokes’s theorem and the 
divergence theorem to show clearly that magnetostatic fields have zero 
divergence. 

[5] 

(d) Carefully explain the consequences of the relationship which you de¬ 
rived in part (c). 
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1996 Test2 Question4 (15 MARKS) 


The diagram below approximates the hysteretic loop for a 100cm 3 sample 
of magnetic material that is the core of a multitap transformer operating at 
50Hz and converting 3A (rms) at 220V (rms) to various output voltages. 


(a) Identify 


-400A/m 


[0.8T 

250A/m 


B 

I.2T 


-250A/m 


400AAn 


-0.8T 

-:JL2T 


(i) the coercive force 

(ii) the remanence 


of the transformer core. 


[ 1 , 1 ] 

(b) Starting with the expression for stored magnetic energy, 

dU H = / H • dB dV 
Jv 

explain how and why it is possible to estimate the energy dissipated in 
a unit volume of a non-linear , isotropic , hysteretic material, given a plot 
of the BH response of the material over a complete cycle of periodic 
magnetic field variation. 

[5] 

(c) Estimate the power loss in watts due to hysteresis in the transformer 
core, and hence state the efficiency of the transformer, assuming that 
2% of the input power is dissipated in the copper windings and that 
the only other loss is due to hysteresis. 


[5,3] 














Specimen Exam Questions 


Q1 (20 MARKS) 

(a) State Faraday’s Law of Electromagnetic Induction, and give a reasoned 
explanation, in terms of magnetic flux linkage, of what the Law means. 

[ 8 ] 

(b) A large electromagnet has circular poles, each of radius 25 cm as shown 
in this diagram. You may assume that the magnetic field is uniform between 
these poles in the z-direction as shown, and zero elsewhere. Use Faraday’s 
Law in the vicinity of a loop of radius r to show that the electric field intensity 
vector between the poles has circumferential component only, and that this 
component is given by E#, where 



«1 C3 


(c) Explain why the E-field has greatest magnitude where r = 25 cm. 

[3] 

(d) It is proposed to raise the magnitude of the B-field from zero to some 
maximum value, B max . This is to take place over a period of 250 ns, in such 
a way that the air between the poles will not break down. At sea level air 
breaks down at about 3000 kV m -1 . Find the maximum value to which the 
B-field can be raised under these circumstances. 

[ 5 ] 
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Q2 (17 MARKS) 

(a) A coaxial cable consists of a solid inner cylinder of radius a and a thin 
outer cylindrical conductor of radius b. The inner cylinder carries a net 
charge of density p Cm -3 . The length of the cable is very much greater 
than the radius of either cylinder. Use Gauss’s Law, with an appro¬ 
priate gaussian surface and a suitable co-ordinate system, to find the 
electric flux density vector both inside and outside the inner cylinder. 
Then calculate the divergence of electric flux density both inside and 
outside the inner cylinder (but inside the outer cylinder), and comment 
on your results. 


[5] 

(b) Current flows in the outer conductor in one direction and back in the 
other direction in the inner conductor. Use Ampere’s Circuital Law to 
explain clearly an advantage of carrying current in a coaxial structure 
rather than, for example, simply using two separate wires. 

[ 2 ] 

(c) The inner cylinder is at a potential of 2V with respect to the earthed 
outer cylinder. The radii of the two cylinders are a = 2.479 mm for the 
inner, and b = 6.738 mm for the outer. Use this information to solve 
Laplace’s equation in the region between the two cylinders, in order to 
arrive at the following formula for the distribution of the potential at 
any radius r: 

$ = —2 In |r| — 10 


[5] 


(d) From the answer to part (c), find the radii of the cylinders of equipo- 
tential at 0.5V, 1.0V and 1.5V. Draw a scale diagram (10 times actual 
size) of the cross section of the conductors and of the equipotentials. 
Use this diagram to find the capacitance per metre of the coaxial cable, 
given that the separating medium has a relative permittivity of e r = 2. 


[ 5 ] 
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Q3 (20 MARKS) 

(a) Explain briefly what is meant by: 

1. Conduction Current 

2. Convection Current 

3. Displacement Current 

In each case you should also supply an equation which defines the stated 
quantity. 

[ 6 ] 


(b) Show how the concept of displacement current is used in the derivation 
of a full statement of Ampere’s Circuital Law: 


f H ■ dl = ■ dS + ^ J s ~D ■ dS 


[7] 


(c) Use the law you derived in part (b) to show that the displacement current 
density within the capacitor in the circuit below is equal to the conduction 
current density along the connecting wires of the circuit. 


Vo SW.(wt> 



[ 7 ] 



SPECIMEN EXAM QUESTIONS 


150 


Q4 (17 MARKS) 

(a) Given that the magnetic field due to a current 1 flowing in a circular 
loop of radius a that lies in the xy-plane is purely z-directed and that, 
at any point z along the axis of the loop, it has value 


* 2 (a 2 + r 2 )§ 

show that the magnetic field of a solenoid of n turns per metre, carrying 
current /, is 


H z = 


n/ 

V(f) 2 + i 


where a is the radius of the solenoid and b is half its length. 


[ 4 ] 


(b) Using the result from part (a), discuss the magnitude and direction of 
the magnetic field of a solenoid of infinite length and of radius a: 

(i) on the axis of the solenoid; 

(ii) within the solenoid but off the axis; 

(iii) outside the solenoid. 

In each case, provide a coherent proof for your stated conclusions and 
explain any assumptions which are made in your analysis. 


[ 2 , 2 , 2 ] 


turn over for part (c) 



SPECIMEN EXAM QUESTIONS 


151 


(c) A coiled current-carrier has radius a = 5 mm and length 6=1 cm. It is 
composed of ./V = 1000 turns of fine copper wire that are wound about 
a former of relative permeability /j, r = 125. The wire carries a current 
of I — 500 mA. 

(i) Derive a general expression in terms oi a, b, N, fi a , fi r and I for 
the energy stored per one-metre length of such a coil, and hence 
find the energy stored in the magnetic field of the particular coil 
described; 

[ 4 ] 

(ii) Show that the inductance of the particular coil is very nearly 1 H. 

[3] 
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Q5 (20 MARKS) 


(a) Use Maxwell’s equations to show that, for an electromagnetic wave trav¬ 
elling in free space in the ^-direction, 


and 


dE x dH y 

dz ^ dt 

dH y dE x 
dz ~ 6 dt 


[5] 


(b) Show that, if E and H both vary sinusoidally, the equations in part (a) 
give rise to the one-dimensional Helmholtz equation: 


cPE x 
dz 2 


—u?[itE x 


Show that a solution of this equation is 


E x = Cl e- Jfc2 + c 2 C kz 


where Cj and c 2 are constants and k = oJy/JH 


[5] 


(c) Hence show that jf- = for the positively-travelling wave solution to 
Maxwell’s equations in free space. 

[5] 


(d) Make a sketch of an electromagnetic wave travelling in the ^-direction in 
free space. Show clearly the directions of E and H at all points along the 
wave, and indicate the ratio of the magnitudes of E and H as well. 


[ 5 ] 
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Q6 (17 MARKS) 

(a) Explain why a high-frequency signal being beamed to a satellite (which 
would typically have a wavelength of a few centimetres) is well-described 
over most of its path as a plane wave in free space. 

[ 2 ] 

(b) Write down Maxwell’s equations for free space, and use them to derive 
the one-dimensional wave equation that governs the x-component of 
electric field, E x , in a plane wave propagating through free space in the 
z-direction. Explain each step in your derivation with care. 

[7] 

(c) Show that 

E x = COS (t - y/fi 0 t 0 z) + COs(f -I- y/Jj^tlz) 
is a solution to the equation that you derived in part (b). 

[ 2 ] 

(d) With reference to the solution suggested in part (c): 

(i) show that a signal beamed to a satellite propagates at a velocity 
close to that of light in a vacuum; 

[3] 

(ii) describe the effects that may be expected when the signal leaves 
the atmosphere of the earth. 


[3] 
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Q7 (20 MARKS) 


(a) Starting from the standpoint of field theory, show that the voltage drop 
across a purely inductive element through which a current i(t) flows is given 
at any instant by 


v = 


L 


di(t) 

dt 


as assumed in circuit theory. Include in your explanation a working definition 
for the quantity L. 


(b) Two coils, through which currents ii(t) and i 2 (t) 
flow respectively, are placed in close proximity to each other 
as shown in this diagram. State clearly what is meant by 
magnetic flux linkage , and mutual inductance , M, and 

show that the voltage induced across the first coil by the 
current in the second is 

Vi 2 = 


[5] 



[5] 



(c) Derive an expression for the inductance per unit length of a coaxial cable. 
Show all your working clearly. You may assume that the coaxial cable consists 
of two excellent conductors of radii a and b, and that it has negligible internal 
inductance. The separating medium has permeability fi. 

[5] 


(d) Using your answer to part (c). show that the external inductance per 
metre of a coaxial transmission line in which the inner and outer conductors 
are separated by a medium of relative permeability /x r =1.005 is 2nH if the 
outer conductor has inner radius 1.01 times the radius of the inner conductor. 


[ 5 ] 


TRTB7R5V 
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Q8 (20 MARKS) 


(a) Either by dimensional analysis, or else arguing from the standpoint of 
field theory, show that the expression 

J = aE 

is a statement-of Ohm’s Law that leads directly to the familiar form in circuit 
analysis, namely 

V = IR 


if a lumped element of length /, conductivity cr, cross-sectional area A and 
resistance 



is considered. 


[4] 


(b) Explain what is meant by continuity, and why it is adequately described 
in a good conductor by 


V J = - 


dp 


dt 


[3] 


(c) Use Ohm’s Law, Maxwell’s equations and the principle of continuity to 
show that a sinusoidally-varying current density in a material of permeability 
p and conductivity a is governed by the equation 

V 2 J = jutpcrj 


[ 6 ] 


(d) Solve the equation in part (c), and hence explain what is meant by the 
concept of skin depth in a good conductor. 


[7] 







































